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ABSTRACT 


The  extensive  cost  to  thoroughly  compare  new  radar  sensor  systems  is  a  problem  in 
today’s  military.  Due  to  the  shrinking  defense  budget,  the  opportunity  to  replace  dated 
sensor  systems,  with  technologically  advanced  systems,  seldom  arises.  Current  funding 
levels  no  longer  support  long  term  evaluations  of  sensor  system  performance.  Die 
development  of  new  methods  to  measure  system  performance  is  crucial  in  determining  the 
best  sensor  system  among  many  alternatives.  Computer  simulation  is  one  method  of 
conducting  additional  trials  to  characterize  sensor  system  performance.  Computer 
simulation  can  aid  decision  makers  in  selecting  the  sensor  system  that  best  meets  the  needs 
of  the  current  military  force  structure.  The  cost  of  simulation  modeling  is  considerably  less 
than  repeated  testing  of  the  real  sensor  system.  This  research  investigates  the  feasibility  of 
developmg  a  computer  simulation  of  a  radar  sensor  system.  The  scope  of  the  research 
includes  computer  modelmg  of  the  detection  process  and  an  evaluation  of  model  output. 
This  simulation  model  is  an  initial  step  to  emphasize  the  power  of  computer  simulation. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made,  within  the 
time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors,  they 
cannot  be  considered  validated.  Any  application  of  these  programs  without  additional 
verification  is  at  the  risk  of  the  user. 
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EXECUTIVE  SUMMARY 


Faster  and  more  complex  radar  sensor  systems  are  being  proposed  and  developed. 
Determining  the  suitability  of  these  new  radar  sensor  systems  as  replacements  for  cunrent 
systems  requires  extensive  testing  of  the  new  system’s  performance.  Thorough  evaluation  of 
new  radar  sensor  systems  incurs  high  cost  and  an  excessive  amount  of  time.  The  declining 
defense  budget  has  an  impact  on  the  quality  of  testing  performed  on  a  given  system  while 
force  reduction  affects  the  number  of  resources  available  to  conduct  complete  system 
evaluation. 

The  quality  of  a  sensor  system  directly  influences  decisions  made  by  tactical 
commanders.  Sensor  systems  that  provide  timely  and  accurate  information  clarify  a  tactical 
commander's  appraisal  of  events  occurring  in  his  operating  environment.  Propier  assessment 
of  the  operating  environment  reduces  the  chance  of  committing  crucial  errors  that  lead  to 
incidents  such  as  fratricide.  A  sensor  system  that  undergoes  rigorous  examination  before 
entering  operational  service  will  reduce  the  chance  of  an  error  in  judgment  by  a  tactical 
commander  using  the  sensor  system. 

Since  limited  funds  reduce  the  extent  of  testing  new  systems,  new  ways  of  augmenting 
the  traditional  testing  methods  need  exploration.  This  thesis  focuses  on  computer  simulation 
modeling  as  an  affordable  and  practicable  procedure  for  augmenting  testing  and  evaluation 
of  new  radar  systems.  Simulation  modeling  of  a  radar  system  allows  for  unlimited  runs  and 
instant  playback  of  test  scenarios.  The  design  of  every  scenario  emulates  the  actual  test 
environmental  conditions  and  interactions  between  the  radar  system  and  its  targets. 
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This  research  led  to  the  development  of  a  computer  program  that  simulates  a  radar 
sensor  system  attempting  to  detect  a  target.  Using  the  operating  characteristics  of  a  given 
radar  system  as  the  inputs,  the  program  determines  how  quickly  the  simulated  radar  is  able  to 
acquire  an  incoming  target.  The  program  outputs  the  amoimt  of  time  required  for  the 
simulated  radar  to  detect  the  target.  To  compare  multiple  radar  sensor  systems,  the  various 
operating  characteristics  are  successively  fed  into  the  program.  The  radar  system  that  detects 
the  target  the  earliest  is  the  preferred  radar  system.  The  effect  on  radar  system  performance 
due  to  varying  its  operating  characteristics  is  quickly  and  easily  visualized.  A  graphical 
display  of  every  scenario  illustrates  the  orientation  and  progress  of  the  engagement  as  it 
evolves. 

This  research  demonstrates  one  possible  method  of  augmenting  traditional  radar 
system  testing.  As  evidenced  by  the  results,  computer  simulation  is  capable  of  providing 
meaningful  results  to  the  user.  By  accurately  describing  the  functioning  of  a  radar  system  by 
use  of  computer  simulation,  strengths  and  weaknesses  of  the  radar  system  are  highlighted 
before  conducting  actual  radar  system  tests.  Correcting  known  deficiencies,  as  identified  by 
computer  simulation,  is  a  more  appropriate  use  of  limited  funds  and  resources. 
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1.  INTRODUCTION 


Faster  and  more  complex  radar  sensor  systems  are  being  proposed  and  developed. 
These  new  radar  sensor  systems  are  candidates  for  replacing  the  current  inventoiy  of 
aging  radar  sensor  systems.  There  are  two  main  questions  to  address  before  replacing  an 
aging  radar  system.  First,  which  new  radar  sensor  system  best  meets  the  minimum 
performance  criterion?  Second,  is  the  potential  replacement  radar  sensor  system  actually 
better  than  the  existing  radar  sensor  system? 

For  the  tactical  commander  embarked  in  a  Naval  ship,  a  radar  sensor  system  is  the 
main  asset  for  building  situational  awareness  and  making  the  state  of  the  environment 
visible.  The  radar  sensor  system  lets  the  tactical  commander  quickly  locate  and  identify 
formations  of  friendly,  enemy,  and  neutral  shipping.  Locations  of  friendly  aircraft  are 
discernible  as  well  as  threat  aircraft  locations  and  maneuvers.  In  both  the  surface  and  air 
environments,  the  tactical  commander  can  assign  friendly  forces  to  investigate 
unidentified  radar  targets  to  further  clarify  the  state  of  the  environment.  Earlier  and  more 
accurate  visualization  of  the  tactical  commander’s  operating  environment  will  directly 
affect  the  safe  operation  of  fnendly  forces  in  a  hostile  environment. 

The  fundamental  method  of  comparing  alternative  radar  sensor  systems  involves 
conducting  trials  to  ascertain  radar  sensor  system  performance.  Attaining  reliable  test 
results  requires  strict  adherence  to  conditions  in  the  test  environment.  For  example,  to 
reduce  variability  an  acceptable  test  environment  requires  the  same  location  and  weather 
conditions.  The  sea  state  at  the  time  of  testing  is  another  variable  that  will  affect  test 
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results.  The  cost  of  holding  these  three  variables  constant  during  testing  is  prohibitively 
high.  Cost  is  not  only  incurred  from  installing  the  competing  radar  sensor  systems  on 
board  a  test  vessel,  but  in  attaining  the  necessary  assets  to  act  as  “targets”  for  the  radar 
sensor  systems  to  detect  and  track. 

Another  disadvantage  of  the  ftmdamental  test  method  is  the  extensive  amount  of 
time  it  takes  to  evaluate  the  multiple  radar  sensor  systems.  To  test  an  air  search  radar  for 
a  ship  requires  the  test  platform  to  transit  from  the  pier  to  the  designated  operating  area. 
Aircraft  fly  out  from  their  home  base,  which  can  be  minutes  to  hours  away  from  the 
ship’s  operating  area.  After  coordination  of  both  units  for  the  test  scenarios,  and 
assuming  no  failures  of  either  platform  or  the  radar  sensor  system  that  is  being  tested,  up 
to  an  hour  of  time  may  have  elapsed.  If  the  aircraft  has  enough  fuel  to  run  trials  for  an 
additional  hour,  perhaps  four  good  runs  are  achievable  before  returning  to  base  for  fuel 
and  maintenance.  Therefore,  only  a  small  amount  of  time  is  truly  available  to  test  a  radar 
sensor  system,  and  achieving  a  full  range  a  performance  testing  is  virtually  impossible. 
Due  to  the  high  cost  and  limited  test  time  available,  the  fundamental  test  method 
provides  only  a  rough  estimate  of  the  true  performance  of  the  radar  sensor  system. 

An  enhancement  to  the  fundamental  method  of  testing  is  computer  simulation 
trials  of  the  competing  systems.  Computer  simulation  provides  the  following  advantages: 

•  Experiment  control  through  the  use  of  set  environmental  standards. 

•  The  variation  of  test  criteria  in  the  evaluation  of  radar  system  performance. 

•  Comparison  of  multiple  alternative  radar  sensor  systems. 
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Cost  and  time  savings  make  computer  simulation  an  attractive  alternative  to  the 
fimdamental  method  of  radar  system  testing.  [Ref  1,  p.  115] 

The  focus  of  this  thesis  is  to  develop  a  computerized  simulation  of  a  radar’s 
detection  process  as  a  means  to  compare  different  radar  systems’  performance  in  a 
tactical  environment.  Of  primary  interest  are: 

•  The  degree  of sophistication  of  the  model 

•  Structure  of  the  model 

•  Required  inputs  to  the  model 

This  work  is  thus  a  first  step  towards  the  use  of  computer  simulation  to  evaluate  different 
radar  systems’  performance  under  the  same  test  and  operational  conditions. 

This  thesis  starts  off  with  a  discussion  of  sensor  systems.  Next,  is  a  discussion  of 
the  simulation  development  in  Chapter  .HI,  followed  by  a  mock  test  scenario  in  Chapter 
IV.  Results  of  the  mock  test  are  presented  in  Chapter  V.  Conclusions  and 
recommendations  conclude  this  study  in  Chapter  VI. 
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IL  SENSOR  SYSTEMS 


A.  BASIC  SENSOR  SYSTEMS 

A  sensor  is  a  device  that  responds  to  actions  occurring  within  its  vicinity.  A 
sensor  system  consists  of  a  response  processor  that  is  connected  to  a  sensor.  The  sensor 
sends  signals  to  the  response  processor  which  vvill  accumulate  data  to  interpret  what  is 
happening  in  the  vicinity  of  the  sensor.  There  are  many  different  types  of  sensor  systems 
that  may  be  applied  in  a  variety  of  ways  in  the  commercial  and  military  world.  They 
range  from  simple  systems,  such  as  photosensitive  cells  that  will  turn  on  a  light  when 
darkness  arrives,  to  a  complex  fire  control  radar  that  will  track  an  inbound  target.  This 
thesis  will  focus  on  complex  radar  sensor  systems. 

Electromagnetic  sensors  may  be  divided  into  three  main  types:  electro-optical, 
laser,  and  radar.  All  sensor  systems  are  designed  to  perform  the  same  task,  albeit  by 
different  methods.  Their  purpose  is  to  aid  the  user  in  determining  the  state  of  the 
environment  by  detecting,  classifying  and  tracking  targets.  A  target  is  a  friendly, 
unknown  or  enemy  unit.  Targets  take  the  appearance  of  personnel  for  an  electro-optical 
sensor  system,  armored  vehicles  for  a  laser  sensor  system  or  aircraft  for  a  radar  system. 
[Ref  2,p.  1] 

Regardless  of  the  method  for  evaluating  the  state  of  the  environment,  the  same 
basic  process  will  be  followed  for  all  sensor  systems  from  initial  search  for  targets  to 
target  destruction.  The  first  event  in  this  process  is  the  detection  of  targets. 
Accomplishing  target  detection  involves  searching  the  environment  using  the  sensor 
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system  xmder  manual  or  automatic  control.  Surveillance  effort  during  the  detection  phase 
focuses  on  all  target  types  in  the  environment. 

Following  the  detection  event  comes  the  proper  identification  of  targets  in  the 
environment.  The  objective  is  to  correctly  identify  and  classify  the  target  as  a  friendly, 
unknown  or  enemy.  Accurate  identification  of  targets  reduces  the  risk  of  fratricide. 
Once  a  target  has  been  identified  it  then  is  placed  into  a  specific  class  of  weapon  system. 
Classification  of  the  weapon  system  allows  the  tactical  commander  to  prepare  a  course  of 
action  in  the  event  that  the  target  turns  hostile. 

If  the  target  is  one  of  interest,  the  next  event  for  the  sensor  system  is  acquisition. 
Difficulties  in  acquisition  of  an  elusive  target  may  arise  due  to  target  maneuvering  to 
avoid  detection  or  target  employment  of  deception  techniques.  In  this  case,  flexibility  in 
changing  sensor  system  operating  parameters  will  increase  the  chance  of  acquiring  the 
target  of  interest. 

The  process  concludes  with  the  target  tracking  event.  After  successful  acquisition 
the  sensor  system  will  track  the  target  using  a  recursive  routine.  Information  related  to 
the  target  such  as  its  speed  and  location  receives  continuous  updates.  An  accurate 
tracking  routine  predicts  with  precision  the  future  location  of  the  target  from  the  targets 
past  dynamics.  Figure  2.1  depicts  a  simple  flow  diagram  of  events  that  every  sensor 
system  will  follow. 


Figure  2.1  Event  processing  of  the  sensor  system 
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B.  DEFINITION  OF  TERMS 

Common  terminology  describes  most  elements  of  sensor  system  functionality. 
Although  the  focus  of  this  thesis  is  radar  systems,  with  minor  changes  the  terminology  is 
applicable  to  other  sensor  systems.  Figure  2.2  shows  a  block  diagram  of  a  simple  radar 
system.  The  block  diagram  illustrates  the  basic  path  a  received  signal  will  follow  within 
the  radar  system. 


Figure  2.2  Radar  sensor  system  block  diagram.  After  (Ref.  3]  [Ref.  2,  p.  82]. 

During  radar  system  operation  the  system  transmits  and  receives  energy  pulses. 
The  pulse  duration  t  is  the  length  of  time  that  power  is  transmitted  during  a  single  pulse. 
The  peak  transmitted  power  level  Pp,  pulse  repetition  frequency  PRF  of  the  system  and  t 
determine  the  average  power  transmitted,  given  by  Equation  2.1. 

P^^P^rPRF  (2.1) 

Figure  2.3  depicts  the  calculation  of  the  range  from  the  sensor  to  the  target.  Time  t 

represents  the  amount  of  time  it  takes  for  the  emission  and  return  of  the  transmitted 
energy  pulse.  Since  the  pulse  propagates  through  the  atmosphere  at  the  speed  of  light  c,  a 
simplistic  calculation  of  range  is  given  by  Equation  2.2.  [Ref  2,  p.  93]  [Ref  2,  p.  56] 

The  presence  of  noise  in  the  radar  receiver  degrades  the  chance  of  the  transmitted 
energy  pulse,  or  signal,  being  classified  as  a  detection.  Noise  comes  mainly  from  three 
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sources:  thermal  noise  associated  with  the  electronic  circuits  in  the  radar’s  receiver,  noise 
from  the  first  stage  in  the  radar  receiver,  and  atmospheric  noise  detected  by  the  radar’s 


Transmitted  Pulse  Returned  Pulse 


time  t 

Figure  2.3  Transmitted  and  returned  pulses. 

From  |Ref.  2,  p.  57] 

R  =  j  (2.2) 

receiving  antenna.  Additional  sources  of  noise  come  from  weather  conditions,  terrain, 
and  other  radio  frequency  emitters.  Upon  entering  the  receiver,  the  received  signal  is 
routed  to  an  intermediate  frequency  amplifier  to  increase  the  chance  of  a  weak  received 
signal  being  classified  as  a  detection.  The  amplification  of  the  return  signal  also  causes 
the  amplification  of  the  noise  signal.  Although  strengthening  of  the  return  signal  through 
amplification  is  helpful  for  signal  processing,  amplification  does  not  ensure  that  the 
minimum  detectable  return  signal  will  be  separable  from  the  noise  inside  the  radar  sensor 
system.  [Ref  4,  p.  3] 

When  the  receiver  output  signal  exceeds  a  predetermined  threshold,  a  detection 
has  occunred.  The  threshold  or  bias  level  can  be  exceeded  by  a  signal  consisting  of  noise 
alone.  The  higher  the  bias  level  setting,  the  less  likely  noise  alone  will  exceed  the  bias 
level  and  trigger  a  detection  Avithin  the  radar  system.  If  noise  alone  causes  a  detection 
then  the  false  detection  is  referred  to  as  a  false  alarm.  [Ref  4,  p.  21]  [Ref  5,  p.  71] 


8 


Gain  is  a  significant  parameter  indicative  of  radar  sensor  system  performance. 
Gain  is  an  effectiveness  measure  regarding  the  ability  of  the  radar  antenna  to  focus  its 
transmitted  energy  in  a  specified  direction.  Antenna  gain  G  as  a  function  of  the  radar 
systems  wavelength  X,  and  aperture  area  A  is  shown  in  Equation  2.3  .  [Ref  2,  p.  3] 


The  aperture  area  is  the  area  of  the  radar  that  receives  the  returned  power  density  that  has 
been  reflected  from  the  target  [Ref  6,  p.  7].  Equation  2.4  gives  the  aperture  as  a  function 


of  the  antenna  gain  and  wavelength. 


Pulse  integration  involves  summing  the  voltages  registered  by  incoming  signals  to 


determine  the  occurrence  of  target  detection.  The  detector  does  not  discriminate 


between  signal  plus  noise  voltages  and  noise-only  voltages.  Therefore,  it  is  possible  to 
sum  only  noise  pulses  to  compare  against  a  threshold  for  determination  of  a  target 
detection.  A  signal-to-fluctuating  noise  ratio,  known  as  standard  deviation  of  noise  or 


root  mean  square  (RMS)  noise,  is  compared  to  a  threshold  value.  If  the  signal-to- 
fluctuating  noise  ratio  exceeds  the  threshold  value  then  a  target  has  been  detected.  The 
advantage  to  pulse  integration  is  that  the  summed  pulses  provide  a  better  indication  of 
the  presence  of  a  target  since  noise  oscillates  aroimd  a  mean  value. 


If  noise  maintained  a  constant  value  over  all  returned  signals,  then  target 
detection  would  consist  simply  of  subtracting  the  noise  level  from  the  returned  signal  and 
determining  if  this  new  signal  is  greater  than  the  threshold.  If  the  signal  minus  constant 
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noise  is  greater  than  the  threshold,  then  a  detection  is  registered.  Therefore,  on  any  given 
pulse  an  accurate  assessment  of  the  presence  of  a  target  is  possible  and  pulse  integration 
is  not  necessaiy. 

In  the  case  of  fluctuating  noise,  the  noise  cannot  be  as  easily  extracted  from  the 
total  signal  since  the  magnitude  of  noise  varies  on  each  pulse.  Accepting  a  signal-to- 
noise  ratio  greater  than  the  threshold  in  this  situation  constitutes  a  false  detection. 
Possibly,  a  large  portion  of  the  total  signal  consists  of  noise  since  noise  can  add  enough 
strength  to  a  signal  to  exceed  the  threshold  setting.  By  integrating  pulses,  the  true  signal 
could  be  averaged  over  n  pulses  to  attain  a  positive  target  detection.  As  a  consequence, 
increasing  the  number  of  pulses  integrated  increases  the  chance  of  positively  detecting  a 
target  at  the  expense  of  a  lengthening  time  until  detection.  [Ref  3,  p.  25]  [Ref  5,  p.  73] 

Transmitted  energy  from  a  radar  sensor  system  experiences  many  types  of  signal 
attenuation,  or  loss,  from  the  time  it  leaves  the  antenna,  travels  through  the  atmosphere, 
and  ultimately  returns  to  the  receiver.  Some  examples  are  scan  distribution  loss,  target 
fluctuation  loss,  integration  loss  and  atmospheric  loss.  System  loss  is  an  aggregate  of 
loss  due  to  radar  system  hardware.  Atmospheric  loss  results  from  absorption  of  signal 
energy  by  the  atmosphere  during  signal  propagation.  The  effect  of  loss  on  radar 
performance  is  the  reduction  in  signal-to-noise  ratio.  Reducing  the  signal  to  noise  ratio 
decreases  the  chance  of  detecting  a  target.  [Ref  4,  pp.  15-16] 

C.  RADAR  EQUATION 

The  radar  equation  forms  the  basis  for  determination  of  target  range.  The 
maximum  range  of  a  radar  sensor  system  is  a  function  of  three  fundamental  parameters; 
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transmitted  power,  antenna  gain,  and  receiver  sensitivity.  Equation  2.5  defines  the 
functional  form  of  the  radar  equation  with  the  signal-to-noise  ratio  as  a  function  of  radar 
system  and  environmental  parameters.  From  Equation  2.5  it  is  obvious  that  the  signal-to- 
noise  ratio  is  inversely  proportional  to  target  range  raised  to  the  fourth  power. 


where 

Pave  =  average  transmitted  power 
A  —  anteima  aperture 
G  =  antenna  gain 

target’s  radar  cross  section 
E  =  integration  efficiency 
n  =  number  of  pulses  integrated 
R  =  range  to  target 

^  noise  ~  TOOt  mean  squared  thermal  noise  power 
==  root  mean  squared  clutter  power 
Lsys  =  system  losses 
Latm  =  atmospheric  losses 

Antenna  aperture  is  the  area  of  the  anteima  that  receives  the  focused  returned  energy 
pulse.  Integration  efficiency  measures  how  well  the  pulse  integrator  operates,  where  a 
flawless  integrator  has  an  efficiency  of  one.  Clutter  causes  distortion  of  the  reflected 
signal  due  to  excess  scattering  of  the  signal  energy  or  absorption  by  some  other  medium. 
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Equation  2.6  demonstrates  that  the  signal-to-noise  ratio  at  time  of  detection  determines 
target  range.  [Ref.  2,  pp.  89, 95]  [Ref  4,  p.  28] 


R  = 


RCS^ 


r  ~ 

\6k^ 

1 _ 

(2.6) 


D.  PARAMETERS 

The  radar  equation  in  Equation  2.6  illustrates  that  average  power,  gain  and 
aperture  define  the  basis  of  radar  sensor  system  operation.  Although  these  parameters 
are  not  the  only  parameters  that  describe  the  functioning  of  a  radar  sensor  system,  they 
provide  sufficient  information  to  evaluate  the  performance  characteristics  of  the  radar 
sensor  system.  The  list  of  sufficient  parameters  that  determine  radar  range  are  provided 
in  Table  1. 

Calculations  based  on  range  are  a  function  of  the  target  that  the  radar  sensor 
system  attempts  to  detect,  the  performance  characteristics  of  the  radar  sensor  system  and 
atmospheric  conditions.  The  radar  cross  section  refers  to  the  area  of  the  target  that 
reflects  energy  back  to  the  radar  sensor  system  receiver.  The  radar  cross  section  of  a 
target  varies  dependent  on  the  aspect  angle  of  the  target  in  relation  to  the  beam  of  the 
radar.  As  the  target  moves  through  the  atmosphere  its  radar  cross  section  fluctuates  in 
intensity.  Fluctuation  of  the  radar  cross  section  occurs  between  radar  sensor  system 
interrogations  of  the  target.  Target  interrogation  happens  on  a  pulse-to-pulse  or  scan-to- 
scan  basis.  [Ref  4,  pp.  37-38, 60] 
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ft  1  DHloniHi 

UNITS  1 

Radar  Cross  Section 

_  CJpr'c 

Gain 

G 

dB 

Wavelenath 

X 

meters 

Power 

P 

.  n\}/f  . . 

Watts 

Anerture 

A 

Table  2.1  Defining  parameters  for  radar  sensor  system. 

After  [Ref.  2,  p.  25] 

A  target’s  radar  cross  section  is  generally  modeled  as  a  composite  of  a  finite 
number  of  points  that  reflect  transmitted  radar  energy  back  to  the  radar  sensor  system 
receiver.  The  composite  radar  cross  section,  equal  to  the  sum  of  the  individual  cross 
sections,  approximately  follows  the  Rayleigh  distribution  [Ref.  4,  p.61].  Another  case  of 
fluctuating  radar  cross  section  combines  a  dominate  reflecting  point  and  many  lesser 
reflecting  points.  This  one-dominant  scatterer  case  of  target  reflectivity  results  in  the  Chi 
Square  distribution  with  four  degrees  of  freedom.  The  multiple  scatterer  and  one 
dominant  scatterer  form  the  two  main  categories  of  fluctuating  radar  cross  sections.  [Ref. 
2,p.  117]  [Ref  4,  pp.  59-61] 
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ni.  SIMULATION  DEVELOPMENT 


A.  DETERMINISTIC  VS.  STOCHASTIC 

The  cookie  cutter  is  the  basic  detection  model  for  a  radar  sensor  system  as  shown 
in  Figure  3. 1 .  The  cookie  cutter  model  easily  integrates  into  any  radar  sensor  system 
scenario.  The  sensor  has  a  detection  range  of  radius  R  that  covers  the  radar  sensor 
system  for  360  degrees.  In  the  cookie  cutter  model  the  radar  sensor  system’s  location  is 
at  the  center  of  the  circle.  A  detection  occurs  when  the  distance  between  the  sensor  and 
the  target  is  less  than  or  equal  to  the  radius  R.  In  its  simplest  form,  cookie  cutter 
detection  guarantees  a  100  percent  chance  of  target  detection.  Since  the  signal-to-noise 
ratio  fluctuates  from  scan  to  scan,  the  chance  of  detection  at  every  point  within  the 


Detection  range  R 


Figure  3.1  Cookie  cutter  detection  model. 

detection  radius  of  the  radar  sensor  system  cannot  be  100  percent.  A  comparison  of  radar 
sensor  system  performance  using  this  simple  detection  model  will  not  be  sufficient  to 
adequately  evaluate  the  performance  of  radar  sensor  systems.  If  the  range  of  detection  is 
the  performance  measure,  then  the  sensor  with  the  farthest  detection  range  will  always 
detect  the  target  first.  Therefore,  the  simple  cookie  cutter  model  does  not  describe  radar 
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sensor  system  performance  at  the  level  of  detail  needed  for  the  comparison  of 
alternatives.  If  the  simple  cookie  cutter  model  is  expanded  to  incorporate  other  events 
that  affect  the  detection  process,  it  is  possible  to  develop  a  model  for  comparing 
alternative  sensor  systems. 

A  stochastic  detection  model  more  accurately  reflects  changes  that  occur  in  the 
physical  operating  environment  of  the  radar  sensor  system  between  successive  scans. 
These  changes  include,  but  are  not  limited  to,  fluctuation  of  the  target’s  radar  cross 
section,  noise  levels  in  the  radar  sensor  system  receiver  and  noise  generating  atmospheric 
conditions.  The  consequence  of  changes  in  the  physical  environment,  between  radar 
scans,  is  the  oscillation  of  the  strength  of  the  signal  and  noise  at  the  receiver  output. 
Random  events  in  the  environment  affect  the  probability  that  the  radar  sensor  system 
detects  a  target  at  some  range  r<R. 

Evaluation  of  the  stochastic  detection  model  requires  use  the  proper  analytical 
tool.  When  combined,  the  radar  sensor  system  events  described  in  Chapter  2  may  be 
viewed  as  a  process,  and  the  detection  event  as  a  sub-process,  in  the  operation  of  the 
radar  sensor  system.  Simulation  modeling  is  a  valuable  tool  for  process  analysis.  For  a 
radar  sensor  system,  simulation  modeling  allows 

•  Time  based  analysis  of  the  detection  process. 

•  The  ability  to  replay  detection  scenarios. 

•  The  capability  to  explicitly  model  randomness 

•  Analysis  of  factors  that  directly  and  indirectly  influence  the  detection  process 
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B.  PROBABILITY  OF  DETECTION 


When  pulse  integration  occurs  after  the  square  law  detector  it  is  called  post 
detection  integration  and  is  depicted  in  Figure  3.2.  Noncoherent  pulse  integration 


Figure  3.2  S^nal  processing.  After  [Ref.  7,  p.  138] 

happens  when  the  phase  relationship  of  the  received  signal  as  compared  to  the  originally 
transmitted  signal  is  destroyed  in  the  intermediate  frequency  filter.  The  efficiency  E  of 
noncoherent  pulse  integration  is  a  positive  value  less  than  one.  Summation  of  the  pulses 
during  noncoherent  integration  provides  input  to  the  threshold  bias  decision  circuits.  A 
detection  takes  place  if  the  input  exceeds  the  threshold  bias.  [Ref.  2,  pp.  81-85]  [Ref  7, 
pp.  137-139] 

Single-hit  probability  of  detection  is  the  chance  of  detecting  a  target.  An  energy 
signal  transmitted  by  a  radar  sensor  system  contains  at  least  one  pulse.  This  pulse  or 
group  of  pulses  has  a  finite  amount  of  time  to  reach  the  target  and  return  to  the 
receiver  while  the  radar  sensor  system  scans  the  environment.  For  Tj,  the  time  available 
for  detection,  the  relationship  holds  that  »  \l{Pulse  Repetition  Frequency).  The 
radar  sensor  system  receives  a  fixed  number  of  pulses  during  the  time  interval  Tj  and 
pulse  integration  then  takes  place.  Equation  3. 1  defines  the  relationship  between  the 
number  of  pulses  available  for  integration  N,  PRF,  width  of  the  radar  beam  9bw  and 
radar  scan  rate  y.  [Ref  7,  p.  144-145]  [Ref  4,  p.  40] 


y 


(3.1) 
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Single-hit  probability  of  detection  is  a  ftmction  of  false  alarm  time.  The  false 
alarm  time  is  the  amount  of  time  that  elapses  between  instances  of  noise  alone  exceeding 
the  threshold  bias  level.  Stated  another  way,  it  is  the  time  required  for  noise  alone  to 
exceed  the  bias  level  with  probability  of  0.5.  The  probability  of  false  alarm  is  given  in 
Equation  3.2  as  a  function  of  threshold  bias  level  Tj  and  N  ,  the  number  of  pulses 
integrated.  I,  in  Equation  3.2,  is  the  incomplete  gamma  function. 

=  (3-2) 


Equation  3.2  is  closely  approximated  by  Equation  3.3  below. 

p  „  AT/'expC-i;) 


(3.3) 


Knowing  the  probability  of  false  alarm  and  the  number  of  pulses  integrated,  a  threshold 
bias  level  can  be  found.  [Ref  5,  pp.  71,77]  [Ref  4,  pp.  21-23] 

The  Swerling  models  are  the  most  commonly  used  models  to  describe  the  single¬ 
hit  probability  of  detection  for  objects  with  fluctuating  radar  cross  sections.  Swerling 
Cases  1  and  2  describe  targets  whose  radar  cross  section  is  a  composite  of  a  finite 
number  of  scattering  points.  Case  1  specifically  addresses  radar  cross  sections  that 
fluctuate  from  scan-to-scan.  During  the  radar  cross  section  fluctuations  among 
individual  pulses  are  highly  correlated.  Little  to  no  correlation  of  radar  cross  section 
fluctuations  is  evident  during  the  time  between  radar  sensor  system  scans.  This  leads  to 
independent  radar  cross  section  fluctuations  for  the  Swerling  Case  1  model.  Equation  3.4 
gives  the  probability  density  function  for  the  signal-to-noise  ratio  for  a  target  with  a 
fluctuating  radar  cross  section. 
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(3.4) 


X  -  observed  signa!-to-noise  ratio 
x  =  mean  signal-to-noise  ratio 

The  Swerling  Case  2  model  adheres  to  the  same  concept  of  multiple  independent 
scattering  points,  but  considers  pulse-to-pulse  fluctuations  of  the  target’s  radar  cross 
section.  The  Swerling  Cases  3  and  4  models  address  target  radar  cross  sections  with  one- 
dominate  scattering  area.  [Ref  4,  pp.  59-63]  [Ref  8,  pp.  122-126] 

The  maximum  detectable  range  R^,  of  a  target  deflned  by  its  radar  cross  section 
Ones  is  found  from  Equation  2.6  by  setting  the  signal  equal  to  noise,  resulting  in  a  signal- 
to-noise  ratio  of  OdB.  Rg  is  assumed  to  be  the  farthest  range  at  which  the  radar  sensor 
system  could  detect  the  target.  Equation  3.5  defines  the  relationship  of  some  range  R  < 
and  the  detection  signal-to-noise  ratio  x.  [Ref  7,  p.  137] 


(3.5) 


To  express  the  single-hit  probability  of  detection  Pj  as  a  function  of  the  threshold 
bias  Yb,  the  mean  signal-to-noise  ratio  J,  and  the  number  of  pulses  integrated  N  for  a 
Swerling  Case  1  target,  the  exact  formulas  are  given  in  Equations  3.6  and  3.7  for  iV=  1 
and  N  >  1  respectively.  Due  to  the  computational  complexity  of  the  exact  single  hit 
probability  of  detection  in  Equation  3.7,  an  approximation  is  provided  in  Equation  3.8 
[Ref  8,  p.  127]. 
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D  -Yt 

Pd  =  exp  — ^ 

1  + jr 


(3.6) 


When  Nx>  1,  Equation  3.8  is  simplified  to  Equation  3.9where  Pj  is  a  function  of  x,  N 
and  Yf,. 


Also,  Equation  3.5  and  3.9  yield  Equation  3.10  as  a  useable  formula  to  find  the  single-hit 
probability  of  detection.  [Ref.  8,  pp.  124-129] 

Pd  =  (3.10) 

where 
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C.  SIMULATION  DESCRIPTION 

The  simulation  takes  radar  sensor  system  parameters  as  inputs  and  determines  a 
random  range  of  detection  based  on  radar  performance  parameters  and  target  radar  cross 
section.  The  basis  for  building  the  simulation  model  is  the  Naval  Postgraduate  School 
Platform  Foundation.  The  NFS  Platform  Foundation  provides  off-the-shelf  building 
blocks  for  simulation  experiments  and  includes  a  graphical  display  of  the  simulation 
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experiment  scenario.  MODSM  II  is  the  object  oriented  language  used  to  program  the 
modules  of  the  NPS  Platform  Foundation.  The  structure  of  the  NPS  Platform  Foimdation 
consists  of  an  imaginary  line  that  divides  the  generic  modules  common  to  all  simulations 
(below  the  line)  and  scenario  specific  modules  designed  by  the  user  (above  the  line). 
Above-the-line  modules  utilize  the  basic  functions  of  below-the-line  modules.  Appendix 
A  includes  a  further  discussion  of  the  NPS  Platform  Foundation.  [Ref  9,  pp.  3-5] 

The  simulation  model  focuses  on  the  detection  process  of  the  radar  sensor  system, 
although  a  complete  scenario  from  detection  to  target  destruction  can  be  run.  The  model 
provides  statistics  concerning  the  time  until  target  of  detection  and  about  the  target’s 
radar  cross  section.  An  example  of  a  scenario  consists  of  a  ship  with  a  radar  sensor 
system  which  is  detecting  an  inbound  aircraft  flying  through  the  radar  sensor  system’s 
area  of  coverage  as  illustrated  in  Figure  3.3.  The  farthest  range  ring  from  the  ship  in 
Figure  3.3  is  the  theoretical  maximum  range  of  the  radar  sensor  system.  The  next  range 
ring  indicates  the  maximum  detectable  range  of  the  target  Ro,  based  on  the  target’s  radar 
cross  section  and  the  radar’s  parameters.  The  range  ring  closest  to  the  ship  is  the  random 
detection  range  of  the  target.  This  example  shows  the  case  of  one  target  and  one  radar 
sensor  system.  The  simulation  model  can  evaluate  multiple  radar  sensor  systems  and 
multiple  targets. 

1.  Model  Assumptions 

The  level  of  detail  chosen  for  the  simulation  model  was  designed  to  reflect  the 
minimum  number  of  parameters  needed  to  describe  the  basic  functioning  of  the  radar 
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sensor  system  in  a  generic  environment.  Conditions  that  might  readily  be  observed  in  the 
physical  environment  need  not  be  explicitly  developed  in  the  simulation  model. 


Figure  3.3  Typical  engagement  scenario. 

As  stated  previously,  the  single  hit  probability  of  detection  is  a  function  of  the 

returned  signal.  The  returned  signal  fluctuates  between  successive  radar  sensor  system 

scans.  Signal  correlation  exists  at  the  output  of  the  receiver,  but  only  for  a  short  time 

duration  that  is  approximately  equal  to  the  inverse  of  the  filter  bandwidth  [Ref  4,  p.  39]. 

This  time  corresponds  to  the  amount  of  time  the  target  is  illuminated  during  a  single 

scan.  Decorrelation  of  returned  signals,  from  consecutive  scans  of  the  radar,  arise 

because  the  elapsed  time  between  scans  is  much  greater  than  the  amount  of  time  the 

signals  remain  correlated.  The  resulting  decorrelation  allows  each  scan  to  be  considered 

a  discrete  independent  look  at  a  target.  Since  this  model  focuses  on  the  initial  detection 

of  the  target,  not  repeated  detections,  signal  independence  simplifies  and  expedites 
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mathematical  calculations.  The  penalties  for  assuming  signal  correlation  between  scans 
are  a  need  for  excessive  computer  memory  and  longer  simulation  run  time. 

The  distance  from  the  target’s  initial  location  to  the  sensor  is  set  at  some  value 
greater  than  Rq  and  the  target  always  proceeds  towards  the  sensor.  The  occurrence  of  the 
initial  detection  is  crucial  for  modeling  the  detection  phase  of  the  radar  since  it  is  the 
mechanism  that  determines  the  response  of  the  radar  sensor  system  and  influences  the 
decisions  made  by  the  tactical  commander.  The  probability  of  not  detecting  a  target,  P„d, 
on  a  single  discrete  scan  is  given  by  Equation  3.11.  Since  the  radar  sensor  system  will 
scan  a  target  more  than  one  time,  the  expression  in  Equation  3.12  shows  that  the 
probability  of  not  detecting  a  target  by  scan  /  is  the  product  of  the  probability  of  not 
detecting  the  target  on  all  previous  scans.  Equation  3.13  shows  the  probability  of 
detecting  the  target  by  scan  /.  [Ref.  4,  p.  39]  [Ref  7,  pp.  58-60, 159] 
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The  initial  location  of  a  target  is  assumed  to  be  outside  of  the  range  R^,  assuming 
that  the  coverage  area  of  the  radar  sensor  system  is  in  the  shape  of  a  hemisphere.  A 
target  approaching  a  radar  sensor  system  originates  from  some  point  A,  normally  a 
considerable  distance  from  B,  as  it  proceeds  to  point  B.  In  this  tactical  engagement 
simulation  model,  point  B  will  always  correspond  to  a  specific  radar  sensor  system’s 
location. 
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2.  Model  Inputs 

Model  inputs  depend  on  the  type  of  analysis  to  be  performed.  The  NFS  Platform 
Foundation  has  its  own  set  of  requisite  model  inputs.  Table  3.1  summarizes  the 
additional  inputs  for  running  this  radar  sensor  system  detection  model.  To  characterize 
the  platform  that  the  radar  sensor  system  desires  to  detect  requires  an  input  of  the  radar 
cross  section  of  the  potential  target.  Describing  the  radar  sensor  system  requires  ten 
additional  inputs.  These  attributes  are:  average  power  aperture  size  A,  integration 
efficiency  E,  number  of  pulses  integrated  n,  wavelength  A,  false  alarm  number  n  ,  filter 

bandwidth  B,  and  receiver  temperature  K° ,  system  loss  L^y,,  and  the  time  delay  between 
successive  scans  of  the  same  point  in  space  Ts.  The  clutter  power  and  atmospheric 
loss  Lsys  explain  the  state  of  the  operating  environment  of  the  radar  sensor  system. 

Of  the  attributes  listed  in  Table  3.1,  the  scan  time  is  the  only  parameter  not 
represented  in  the  radar  equation  given  in  Equation  2.5.  The  scan  time  sets  the  proper 
time  interval  for  the  independent  scans  of  the  target.  The  parameters  B  and  K°  relate  to 
the  thermal  noise  power  in  the  receiver  and  is  shown  in  Equation  3.14  [Ref  6,  p.  133], 
As  evidenced  in  Equation  2.5,  thermal  noise  plays  an  important  role  in  the  radar 
equation. 

D.  IMPLEMENTATION 

The  modules  developed  for  this  thesis  were  specifically  designed  to  evaluate  a 
radar  sensor  system’s  detection  phase.  The  additions  to  the  NPS  Platform  Foundation 
include  both  above-the-line  and  below-the-line  modifications.  Below-the-line  changes 
consisted  of  ensuring  the  correct  parameters  were  passed  to  the  new  modules  such  that 
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the  integrity  of  the  NPS  Platform  Foundation  remained  intact.  Other  changes  were  made 
to  add  the  additional  parameters  needed  to  model  the  radar  sensor  system  as  shown  in 
Table  3.1.  Appendix  B  provides  the  source  code  of  the  below-the-line  changes.  The 
above  the  line  changes  were  made  to  develope  the  essence  of  the  radar  sensor  system 
operating  in  a  tactical  environment.  The  structure  is  built  around  an  environmental 
object  as  shown  in  Figure  3.4.  The  environment  transfers  information  between  the  target 
and  the  sensor.  The  sensor  is  an  independent  entity,  but  is  linked  to  its  host  platform 

through  a  virtual  sensor.  Therefore,  referencing  the  sensor  or  the  host  platform  leads  to 
the  same  effect. 


Name 

Symbol 

Units 

Radar  Cross  Section 

^RCS 

Meters^ 

Average  Power 

P 

a\;e 

Watts 

Integration  Efficiency 

E 

Pulses  Integrated 

n 

— 

Wavelength 

X 

Meters 

False  Alarm  Number 

n 

— 

Filter  Bandwidth 

B 

Hz 

Receiver  Temperature 

T 

Scan  Time 

Ts 

Seconds 

Clutter  Power 

—2 

Watts 

System  Loss 

Lsys 

— 

Atmospheric  Loss 

Latm 

— 

Table  3.1  Radar  attribute  inputs  for  the  simulation  model. 

<^Lse^kTB  (3.14) 

where 

yt  =  138  X 10“^ 
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Figure  3.4  Two-way  information  flow  between  environment,  sensor,  and  threat. 

1.  Environment 

The  physical  atmosphere  through  which  the  signal  travels  is  simulated  by  a 
variable  called  the  environment  as  found  in  Appendix  B.  The  environment  is  the  source 
for  determining  when  and  where  a  detection  takes  place.  The  environment  knows  the 
performance  characteristics  of  the  radar  sensor  system,  speed  and  three  dimensional 
location  coordinates  of  the  platform  and  target,  and  the  target’s  radar  cross  section.  The 
role  of  the  environment  is  to  randomly  select  a  detection  time  from  the  cumulative 
distribution  function  (CDF)  for  the  probability  of  detecting  the  target.  Construction  of 
the  CDF  is  based  on  Equation  3.13,  where  each  scan  i  is  determined  by  the  radar  sensor 
system  scan  time  7^.  Every  Tg  seconds  the  radar  will  have  a  discrete  and  probabilistically 
independent  look  at  the  target  with  a  corresponding  probability  of  detection. 

The  process  for  building  the  CDF  requires  determining  the  location  coordinates  of 
Ro  based  on  the  movements  of  the  platform  and  the  target.  Next  the  environment  finds 
the  ranges  and  location  coordinates  of  the  target  for  each  discrete  radar  sensor  system 
scan.  Since  this  thesis  focuses  on  radar  detection,  the  ranges  and  location  coordinates  are 
found  only  while  the  target  travels  inbound  to  the  sensor.  The  cumulative  probability  of 
detection  is  computed  by  an  iterative  process  of  summing  the  probability  of  detection  for 
all  prior  radar  sensor  system  scans.  At  this  point,  for  each  radar  scan,  the  environment 
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has  computed  a  range  where  the  scan  transpires,  the  probability  of  detection  on  the  scan, 
and  the  cumulative  probability  of  detection  on  the  scan. 

The  accumulated  probabilities  form  a  discrete  step  fimction  from  which  a  random 
scan  representing  the  detection  is  drawn.  The  method  for  drawing  the  random  scan 
comes  from  the  inverse  transform  method  of  generating  discrete  random  variates  [Ref  1, 
pp.  469-474].  The  time  associated  with  the  scan  represents  the  additional  time  until 
detection  of  the  target  after  reaching  the  range  Rg.  To  clarify,  represents  the  farthest 
range  that  a  radar  sensor  system,  based  on  its  operating  characteristics,  can  detect  a  target 
with  a  radar  cross  section  of  size  ctrcs-  Calculating  the  time  for  the  target  to  proceed 
form  its  initial  location  until  intercepting  a  hemisphere  around  the  sensor  of  radius  R^  is 
always  a  fixed  time  of  travel.  Therefore,  since  the  target  cannot  be  detected  at  any  range 
greater  than  Ro,  the  only  randomness  comes  in  to  play  after  the  target  reaches  its 
maximum  detectable  range.  The  random  drawn  time  is  added  to  the  time  the  target 
intercepts  Ro  as  determined  by  the  NPS  Platform  Foundation.  When  the  simulation  time 
equals  the  sum  of  the  two  times  a  detection  is  registered.  If  the  graphics  screen  is  utilized 
to  view  the  scenario,  the  user  will  notice  the  range  ring  of  the  sensor  changes  from  green 
to  red  and  the  status  ring  of  the  target  changes  from  yellow  to  blue  at  the  instant  the 
target  is  detected. 

2.  Radar  and  Threat 

The  first  step  for  initializing  the  radar  sensor  system  detection  process  is 
determining  the  range  Ro  from  the  radar  equation.  This  range  is  where  the  signal  to  noise 
ratio  equals  one.  From  Equation  3.5  the  relationship  between  range  and  signal  to  noise 
ratio  has  been  established.  The  radar  object  then  determines  J*,  the  threshold  setting. 
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The  threshold  is  automatically  set  by  a  recursive  routine  Mrhich  can  be  found  in  Appendix 
B  [Ref  4,  p.  486],  At  the  time  of  detection  the  radar  object  gathers  statistics  related  to 
the  signal  to  noise  ratio  that  are  used  for  analysis.  The  threat  object  is  merely  an 
interface  to  provide  information  to  the  radar  via  the  environment  with  out  divulging 
hidden  information  concerning  the  target’s  location,  speed,  and  radar  cross  section. 
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IV.  PROOF  OF  CONCEPT 


The  best  method  of  demonstrating  the  concept  of  using  simulation  as  an  aid  to 
radar  sensor  system  analysis  is  through  a  scenario  emulating  actual  events.  The  scenario 
compares  the  performance  of  two  competing  radar  sensor  systems.  The  expected  time 
until  detection  is  the  measure  of  effectiveness  for  the  scenario.  After  determining  the 
expected  time  of  detection  of  each  system  against  a  target  of  a  given  radar  cross  section, 
the  combined  signal  and  noise  distribution  is  examined.  The  examination  focuses  on  the 
observed  combined  signal  at  the  time  of  target  detection.  From  this  analysis  conclusions 
are  drawn  to  describe  the  functionality  and  compatibility  of  the  radar  sensor  system. 

A.  MOCK  SCENARIO 

The  scenario  involves  the  evaluation  of  the  SPS-30  and  SPS-49  air  surveillance 
radars.  The  SPS-30  and  SPS-49  are  in  competition  to  become  the  replacement  for  the 
current  air  surveillance  radar  system,  which  has  reached  the  end  of  its  service  life.  The 
manufacturers  view  the  chance  to  produce  the  follow-on  to  the  current  air  surveillance 
radar  sensor  system  as  an  extremely  profitable  government  contract  Although  the 
manufacturers  both  claim  to  meet  the  military  specification  for  the  air  surveillance  radar 
sensor  system  design,  only  one  will  be  chosen  as  the  follow  on  to  the  current  air 
surveillance  radar  sensor  system. 

The  high  cost  of  installing  the  competing  radar  sensor  systems  onboard  a  naval 
vessel  for  evaluation  and  limited  availability  of  test  assets  has  led  to  the  use  of  simulation 
modeling  for  initial  exploratory  analysis  of  the  radar  sensor  systems  performance.  After 
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gaining  knowledge  of  the  capabilities  and  weaknesses  of  the  radar  sensor  system  through 
simulation  output  analysis,  resources  can  then  be  allocated  to  investigate  and  test  specific 
areas  of  air  surveillance  system  performance  to  acquire  more  precise  data  for  comparing 
the  two  systems. 

The  scenario  has  the  SPS-30  and  SPS-49  test  being  conducted  in  an  controlled 
ocean  environment.  The  test  involves  the  detection  of  a  target  with  a  radar  cross  section 
size  of  5  square  meters.  The  target  flight  profile  involves  a  constant  altitude  of  500  feet 
and  constant  airspeed  of  500  knots,  commencing  from  a  distance  of  100  nautical  miles 
from  the  ship.  The  metric  for  the  test  is  the  expected  time  of  detection,  calculated  from 
the  time  the  target  has  commenced  its  inbound  run.  The  smaller  the  expected  time  of 
detection,  the  better  the  performance  of  the  air  surveillance  system  under  this 
measurement  criteria.  Operational  characteristics  of  the  two  radar  sensor  systems  is 
fictitiously  provided  by  the  radar  sensor  system  manufacturers.  Table  4. 1  lists  the  input 
parameters  to  the  simulation  model  that  specify  the  radar  sensor  system’s  operating 
characteristics. 


11:1 

SPS-49 

P 

,.5.  MVP _ 

9000 

13000 

A 

8.82 

6.24 

E 

0.98 

0.95 

n 

5 

7 

A _ 

0.15 

0.319 

n’ 

riooooo 

100000 

Ts _ 

6.0 

5.0 

T 

800 

1 B  2000000 

2000000 

Table  4.1  Scenario  inputs.  After  [Ref.  10,  pp.  222, 224] 

The  actual  events  in  a  scenario  are  illustrated  in  Figure  4. 1  through  Figure  4.7. 

When  viewing  the  scenario  on  a  computer  monitor  specific  colors  describe  significant 
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events.  Figure  4.1  displays  the  initial  setup  of  the  ship  and  the  target.  The  small  status 
ring  around  each  platform  indicates  there  has  not  been  a  detection.  Visually  on  a 
computer  monitor  the  initial  color  of  the  status  ring  is  yellow.  In  Figure  4.2,  the  radar’s 
range  ring  representing  the  maximum  detectable  range  of  the  target  is  added  to  the  ship. 
When  viewing  the  scenario  on  a  computer  monitor  the  initial  color  of  the  radar  range  ring 
is  green.  After  the  addition  of  the  radar  range  ring  the  simulation  begins.  The  target 
proceeds  inbound  to  the  ship  in  Figure  4.3,  as  the  ship  proceeds  on  its  course.  Although 
the  target  has  crossed  the  maximum  detectable  range.  Figure  4.4,  its  status  ring  shows 
that  it  has  not  been  detected  by  the  radar  sensor  system.  The  status  ring,  when 
graphically  displayed  on  a  monitor,  retains  its  yellow  color.  Figure  4.5  depicts  the 
detection  event.  Viewed  on  a  monitor,  both  the  status  ring  of  the  target  and  the  radar 
sensor  range  ring  simultaneously  change  color.  Maroon  identifies  the  target  in  a  detected 
status  and  the  ship’s  range  ring  changes  to  red  indicating  the  presence  of  a  target.  As  the 
target  continues  inbound.  Figure  4.6,  its  status  ring  continues  to  display  a  detected  status. 
The  simulation  concludes  when  the  target  overflies  the  ship  as  pictured  in  Figure  4.7. 

B.  SCENARIO  ASSUMPTIONS 

The  modeling  assumptions  for  this  scenario,  and  the  simulation  model  in  general, 
are  devised  to  follow  common  reason  and  decrease  computer  programming  complexity. 
The  first  assumption  for  any  scenario  involving  this  simulation  model  is  that  the  initial 
location  of  the  incoming  target  is  greater  than  Ro.  Also,  target’s  mission  is  to  attack  the 
ship  platform.  Therefore,  the  target  will  proceed  to  the  last  known  coordinate  position  of 
the  ship.  The  target’s  radar  cross  section  fluctuates  in  accordance  with  the  Swerling  Case 
1  target  model.  Provisions  have  been  made  to  easily  add  other  Swerling  models  to 
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Figure  4.1  Initial  setup  of  the  ship  and  the  target. 


32 


Figure  4.2  Sensor  range  ring  established. 


33 


Figure  4.3  Ship  steady  on  course  and  target  inbound. 
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Figure  4.4  Target  crossing  maximum  detectable  range. 
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Figure  4.5  Radar  detects  target. 


Figure  4.6  Detected  target  proceeding  inbound. 
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Figure  4.7  Simulated  engagement  complete. 
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the  current  computer  code.  The  last  general  assumption  stipulates  that  all  range 
calculation  are  given  in  slant  range  to  the  target. 

Assumptions  in  this  scenario  assert  that  the  received  signal  does  not  experience 
the  corrupting  effect  of  backgroimd  clutter.  The  parameters  integration  efficiency,  pulses 
integrated,  false  alarm  number,  receiver  temperature,  and  filter  bandwidth  are  not  true 
operating  values  of  the  two  radar  sensor  systems.  These  values  merely  portray  possible 
actual  values.  The  simulation  run  is  complete  when  the  target  passes  overhead  the  ship. 
Because  this  scenario  evaluates  air  surveillance  radar  sensor  systems  in  a  tactical 
environment,  the  speed  of  the  incoming  target  greater  than  the  speed  of  the  ship. 
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V.  RESULTS 


The  expected  time  imtil  target  detection  was  measured  for  the  SPS-30  and  the 
SPS-40  radar  sensor  systems  based  on  the  input  parameters  from  Table  4.3  The 
procedure  for  collecting  the  data  involved  taking  data  samples  of  the  elapsed  time 
derived  from  the  time  the  target  commenced  its  inbound  run  until  it  was  detected.  The 
samples  from  repeated  runs  of  the  scenario  were  then  combined  to  calculate  an  estimate 
of  the  mean  time  until  detection.  The  point  estimate  is  then  bounded  by  confidence 
intervals  to  further  clarify  the  accuracy  of  the  estimate. 

The  simulation  experiment  is  replicated  until  a  reasonable  estimate  of  the  mean 
time  until  detection  is  found.  The  stopping  rule  for  the  simulation  model  is  the  relative 
error  of  the  mean  time  to  detection.  The  relative  error  7  specifies  the  percent  error  in  the 

estimate  of  the  mean  X  from  the  true  value  of  the  mean  ju.  Equation  5.1  shows  the 
relationship  between  the  mean  values  and  the  relative  error  [Ref  1,  p.  537].  Based  on  the 

’'=V 

magnitude  of  the  relative  error,  the  number  of  replicates  the  simulation  performs  is 
related  to  the  half-width  5  of  the  confidence  interval  for  the  estimate  of  the  mean.  The 
confidence  interval  changes  as  a  function  of  the  number  of  replicates  n.  Each  replicate 
checks  the  size  of  the  confidence  interval  half-width  and  compares  the  half-width  to  the 
relative  error  as  given  in  Equation  5.2.  [Ref  1,  pp.  538-539] 


%)  ^  y 

X{n)\  1  +  7 


(5.2) 
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For  the  mock  scenario  the  relative  error  was  chosen  to  be  1  percent.  Each  radar 
sensor  system  replicate  took  less  than  one  second  to  run.  The  run  time  for  a  particular 
radar  sensor  system  was  less  than  four  minutes.  Based  on  the  mean  time  of  detection  and 
the  given  scenario,  the  SPS-30  performed  better  than  the  SPS-49.  Table  5.1  summarizes 
the  output  results. 


Measure _ 

^S-30 

SPS-49 

Mean  Time  fminutesl 

4.9192 

6.4556 

Confldence  Interval  fCB 

(4.8735.4.9649') 

(6.4121.6.4992') 

Cl  Precision 

0.01 

0.01 

Cl  Width 

0.0913 

rao87i 

Samole  Variance 

0.2882 

0.1317 

Runs 

531 

267 

Table  5.1  Time  until  detection  summary  statistics. 

The  flexibility  of  simulation  modeling  is  further  demonstrated  by  extending  the 

analysis  to  exploring  the  returned  signal  and  noise.  Noise  is  modeled  by  assuming  that  it 

comes  from  the  Rayleigh  distribution  shown  in  Equation  5.3.  The  signal  comes  from  an 

unknown  distribution,  but  is  a  function  of  the  fluctuating  radar  cross  section  of  the  target 

and  the  range  of  the  target  at  the  time  of  detection.  The  distribution  of  signal  plus  noise 

theoretically  comes  from  the  Rice  distribution  given  in  Equation  5.4.  Noise  signals  are 

generated  by  the  inverse  transform  method  for  the  Weibull  distribution  since  the 

Weibull(2,  p)  is  equivalent  to  the  Rayleigh(j3),  where  P  equals  the  square  root  of  two 

times  the  total  noise  power  <jI  [Ref.  1,  pp.  333-334, 490].  [Ref  3] 

-r^ 

P(r)  =  -^exp2"'  (5.3) 

where 

^  ^clutter  ^ noise 
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where 


lo  =  Bessel  function  of  order  0 
A  =  half  wave  of  the  sinusoidal  signal 

The  accumulation  of  signal  plus  noise  data  follows  the  same  method  and  stopping 
rule  as  used  in  collecting  detection  time  data.  The  received  signal  and  fluctuating  noise 
are  computed  then  summed  to  form  the  composite  signal.  Once  the  composite  signal 
data  points  are  gathered,  the  observations  are  compared  to  the  theoretical  distribution  to 
appraise  how  closely  they  reflect  the  theoretical  distribution  of  signal  plus  noise.  Table 
5.2  consolidates  the  results  from  the  two  runs  while  Figure  5.1  plots  a  histogram  of  the 
observed  data  against  the  theoretical  distribution.  All  observed  data  values  have  been 
rescaled  by  a  factor  of  10’’. 


SPS-30 

SPS-49 

Mean  Value 

8.8383 

5.4060 

Confidence  Interval  fCD 

Cl  Precision 

0.035 

0.04 

Cl  Width 

0.6179 

0.4321 

Samnle  Variance 

42.5358 

18.6033 

Runs 

^1712 

1531 

Table  5.2  Signal  plus  noise  summary  statistics. 


2  4 

Signal  +  Noise  Value 
Solicittie  =  Theoretical  Dist 


Figure  5.1  Comparison  of  the  observed  data  to  the  theoretical  distribution. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


TMs  thesis  demonstrates  the  potential  of  simulation  modeling  as  a  tool  to  assist  in 
the  testing  of  radar  sensor  system.  The  results  clearly  state  the  benefit  of  the  ability  to 
replay  a  given  detection  scenario,  as  evidenced  through  the  precision  of  the  mean  value 
estimates  for  the  time  until  detection  and  the  signal  plus  noise  value.  Additionally,  the 
simulation  shows  that  many  replications  of  a  scenario  are  needed  to  determine  a  good 
point  estimate  of  the  true  mean  of  the  variable  under  consideration.  The  fundamental 
method  of  comparing  alternatives  does  not  allow  a  vast  number  of  trials  due  to  the  high 
cost  of  repeating  the  trials.  Fundamental  testing  exercises  very  limited  control  over  the 
environment,  thereby  increasing  variability  in  the  observed  results. 

Simulation  modeling  allows  for  exploratory  analysis  of  sensor  system  processes 
which  can  lead  to  the  discovery  of  new  measures  of  sensor  system  performance.  For 
example,  assume  that  noise  follows  a  logistic  distribution  rather  then  a  Rayleigh 
distribution.  The  effect  on  the  signal  plus  noise  or  signal  to  noise  ratio  may  then  be 
examined.  Converse  to  studying  sensor  systems,  one  might  desire  to  understand  how  a 
new  weapon  system  performs  against  a  particular  enemy  sensor  system.  By  entering  the 
basic  parameters  of  the  enemy  sensor  system,  knowledge  of  the  detectability  of  a 
proposed  friendly  weapon  systems  is  gained. 

This  research  evaluated  a  single  scenario  to  demonstrate  the  feasibility  of 
modeling  radar  sensor  system.  The  ability  exist  to  run  a  multitude  of  scenarios  to 
investigate  possible  weaknesses  in  radar  sensor  system  performance.  Furthering  of 
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research  in  the  area  of  sensor  systems  can  be  done  by  more  accurately  defining  the 
effects  of  the  atmosphere  on  the  propagating  signal,  such  as  ducting  and  multipath. 
Power  intensity  of  a  reflected  signal  could  be  adapted  to  reflect  the  aspect  angle  of  the 
target  relative  to  the  direction  of  the  transmitted  signal.  After  the  desired  level  of  detail 
is  reached,  as  defined  by  the  user,  the  other  processes  of  the  sensor  system  should  be 
modeled. 
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APPENDIX  A.  NPS  PLATFORM  FOUNDATION 


This  appendix  provides  a  brief  overview  of  the  NPS  Platform  Foundation.  The 
Platform  consists  of  modules  designed  to  provide  an  off  the  shelf  generic  simulation 
package  to  run  graphical  military  engagement  scenarios  with  an  emphasis  towards  sensor 
systems.  The  Platform  Foundation  is  written  in  the  object  oriented  language  MODSIM II 
and  the  graphics  are  programmed  in  SIMGRAP H II.  Scenarios  are  built  using  data  input 
files  to  specify  the  types  of  systems  employed,  locations  of  the  systems,  and  system 
performance  characteristics  and  movements.  The  graphical  map  display  depicts  the 
systems  and  their  movements  on  a  100x100  nautical  mile  grid. 

The  basic  components  of  any  scenario  are  the  platforms  and  the  sensors.  The 
platforms  can  be  of  a  specified  type  such  as  a  destroyer,  frigate,  attack  jet,  or  fighter  jet. 
The  platform  is  then  assigned  an  identification  name  that  relates  it  to  a  specific  type  of 
platform,  for  example,  a  platform  type  ‘frigate’  is  created  and  a  platform  is  given  the 
identification  name  of  OHPerry.  The  identification  name  OHPerry  relates  the  platform  to 
a  platform  type  ‘frigate’.  This  building  convention  allows  for  multiple  forces  of  the  same 
type  to  be  present  during  an  engagement,  which  is  comparable  to  real  military  battle 
force  structure.  Sensors  are  an  attribute  of  platform  lypes.  There  is  no  direct  link 
between  a  sensor  and  a  platform  except  through  the  graphical  movements  of  the 
platform.  Whereas  onboard  a  real  naval  vessel  a  fire  control  radar  is  physically  part  of 
the  ship,  in  the  virtual  world  the  only  cormection  between  the  fire  control  radar  and  the 
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ship  is  a  graphic  image  of  the  ship  and  a  range  ring  representing  the  fire  control  radar 
scan  range. 

Sensor  are  linked  to  targets  by  a  virtual  sensor.  Virtual  sensing  is  the  medium 
through  which  information  is  transferred  between  target  and  sensor.  The  virtual  sensor 
maintains  a  queue  of  all  targets  it  is  tracking  and  calculates  the  time  of  detection  based 
upon  movements  of  the  sensor/platform  and  the  target.  Virtual  sensing  keeps  a  platforms 
movements  and  attributes  concealed  from  other  platforms.  Virtual  sensing  is  the  main 
force  driving  any  scenario  involving  sensing. 
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APPENDIX  B.  SIMULATION  CODE 


The  following  code  is  a  combination  of  new  modules  and  critical  additions  to 

existing  NPS  Platfonn  Modules  that  needed  changes  for  running  the  simulation  model. 

All  programming  codes  is  done  in  MODSIU II  while  the  program  is  run  on  a  UNK 

worksmtion  [Ref  1,].  The  first  changes  to  existing  code  was  done  to  Are  modules 

DelKng  and  PLisl.  These  changes  added  the  additional  parameters  needed  to  model  the 

radar  sensor  system  and  target  attributes.  The  change  to  PLrstmod  was  a  single  line 

adding  the  radar  cross  section  as  an  attribute. 

Definition  module  DDetRng.mod: 

FROM  GM  IMPORT  RealmType; 

FROM  ListMod  IMPORT  QueueList; 

FROM  RecIO  IMPORT  RecIOHandlk)bj; 

TYPE 

AttributesRecType  =  RECORD 
AveragePower  :  REAL; 

Aperture  .  REAL; 

EflSciency  ;  REAL; 

PulseInt  :  INTEGER; 

WaveLength  :  REAL; 

FalseAlannNumber ;  REAL; 

ScanTime  ;  REAL; 

ReceiverTemp  ;  REAL; 

FilterBW  :  REAL- 

END  RECORD; 


SensorlnfoRecType  =  RECORD 


SensorTypeName 

Realm 

DefaultRange 

TypeSensor 


:  STRING; 

:  RealmType; 

.REAL; 

:  STRING;  {i.e.  RADAR  OPTICAL,  etc.} 


Attributes 
END  RECORD; 


:  AttributesRecType; 


SensorlOHandleObj  =  OBJECT(RecIOHandleObj[ANYREC- 
END  OBJECT; 


SensorlnfoRecType]) 


49 


SensorMasterListObj  =  OBJECT(QueueList[ANYREC:  SensorInfoRecType]) 

ASK  METHOD  GiveSensor(IN  Name  :  STRING):  SensorInfoRecType; 
END  OBJECT; 

CoupledRangeRecType  =  RECORD 
TgtType  :  STRING; 

SensorType :  STRING; 

Range  :  REAL; 

END  RECORD; 

DetectionRangePairObj  =  OBJECT(QueueList[ANYREC:  CoupledRangeRecType]) 

ASK  METHOD  ReadRangeRecs; 

ASK  METHOD  CoupledRange(IN  TgtType  :  STRING; 

IN  SensorTypeName  :  STRING):  REAL; 

END  OBJECT; 


VAR 

MasterSensorList  :  SensorMasterListObj; 
DetectionRangeOracle :  DetectionRangePairObj; 
END  MODULE. 


Implementation  module  DDetRng.mod:  (Changes  only) 


ASK  METHOD  CoupledRange(IN  TgtTypelD  :  STRING; 

IN  SensorTypeName  :  STRING).  REAL; 


VAR 


CurrentSensor :  SensorInfoRecType; 
BEGIN 


CurrentSensor  :=  ASK  MasterSensorList  TO  CSveSensor(SensorTypeName); 
IF  CurrentSensor.  TypeSensor  =  "RADAR" 

RETURN! ASK  l^dar  TO  FindRo(TgtTypeID,  SensorTypeName)); 

ELSE 


RETURN(-1.0); 
END  IF; 

END  METHOD; 


ASK  METHOD  ReadRangeRecs; 

VAR 

SHRec :  SHRecType; 

SHArray :  SHArrayType; 
i,  j :  INTEGER; 

SensorRec :  SensorInfoRecType; 

CharArray  :  ARRAY  INTEGER  OF  CHAR; 

Rec :  CoupledRangeRecType; 

MasterRangelOHandler :  SensorlOHandleObj; 

BEGIN 

NEW(MasterSensorList); 

NEW(MasterRangeIOHandler); 

ASK  MasterRangelOHandler  TO  ReadRecs("sensor.dat"); 

ASK  MasterRangelOHandler  TO  FindSHRec("SensorTypes",  SHRec); 
FOR  i  :=  1  TO  HIGH(SHRec.OwnedString)  BY  13 
NEW(SensorRec); 

SensorRec. SensorTypeName  :=  SHRec.OwnedString[i]; 
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SensorRec.Realm  ConvertToRealiiiType(SHRec.OwnedString[i+l]); 
SensorRec.DefaultRange  ;=  STRTOREAL(SHRec.OwnedString[i+2]); 

SensorRec.TypeSensor  :=  SHRec.OwnedString[i+3]; 

NEW(SensorRec.  Attributes); 

SensorRec. Attributes. AveragePower  :=  SllRTOREAL(SHRec.OwnedString[i+4]) 
SensorRec.Attributes.Aperture  ;=  STRTOREAL(SHRec.OwnedString[i+5]); 

SensorRec. Attributes-EfBciency  :=  STRTOREAL(SHRec.OwnedStringti+6/); 

SensorRec.  Attributes.PuIseInt  :=  STRTOINT(SHRec.OwnedString[i+7]); 
SensorRec.Attributes.WaveLength  ;=  STRTOREAL(SHRec.OwnedString[i+8]); 
SensorRec.Attributes.FalseAlarmNumber  :=  STRTOREAL(SHRec.OwnedString[i+9])- 
SensorRec.Attributes.ScanTime  :=  STRTOREAL(SHRec.OwnedString[i+10])/60  0- 
SensorRec.Attributes.ReceiverTemp  ;=  STRTOREAL(SHRec.OwnedString[i+l  !])■  ’ 

SensorRec.  Attributes.FaterBW  :=  STRTOREAL(SHRec.OwnedStrinari+121> 

WriteLineC  "); 

WriteLine{"Sensor  "  +  SensorRec. SensorTypeName  +  "  has  realm  "  + 
RealmToStr(SensorRec.Realm)  + "  and  default  range "  + 

REALTOSTR(SensorRec.DefaultRange)+  "  This  sensor  is  a  "+SensorRec.TypeSensor) 
WnteLineC'Attributes  are  Avg  Power:  "  +  REALTOSTR(SensorRec. Attributes.  AveragePower)); 
WriteLine(  Apeture:  "  +  REALTOSTR(SensorRec. Attributes. Aperture)); 

WriteLine("  Efficiency.  "  +  REALTOSTR(SensorRec.  Attributes.Efficiency))- 

WnteLineC  Pulses  Integrated:  "  +  INTTOSTR(SensorRec.  Attributes.PuIseInt));  ’ 

WriteLine(  Wavelength: "  +  REALTOSTR(SensorRec.Attributes.WaveLength))' 

WriteLine("  False  Alarm  Number: " 

+  REALTOSTR(SensorRec.Attributes.FalseAlaitnNumber))- 
WnteLme("  Radar  Scan  Time: "  +  REALTOSTR(SensorRec.Attributes.ScanTime))- 

WriteLineC’  Radar  ReceiverTemp:  ” 

+  REALTOSTR(SensorRec.Attributes.ReceiverTemp)); 
WnteLine("  Radar  FhterBW.  "  +  REALTOSTR(SensorRec.Attributes.FflterBW))- 

ASK  MasterSensorList  TO  Add(SensorRec) 

END  FOR; 

END  METHOD; 

The  modules  D AARON  1. mod  and  I  AARON  1. mod  comprise  the  environment, 
radar,  and  threat  objects.  The  radar  determines  its  threshold  setting  in  this  module  arid  the 


environment  makes  the  link  to  create  the  cumulative  distribution  function  in 

lAARONCDF.mod. 

Definition  module  DAARONl.mod: 

DEFINITION  MODULE  AARONl ; 

FROM  GM  IMPORT  LocationRecType; 

FROM  Plat  IMPORT  PlatformObj; 

FROM  DetRng  IMPORT  SensorInfoRecType; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  MathMod  IMPORT  POWER; 

FROM  SimpleStats  IMPORT  StatObj; 

FROM  RandMod  IMPORT  RandomObj; 
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TYPE 

ThreatObj  =  OBJECT 

ASK  METHOD  Objinit; 

ASK  METHOD  CalculateRCS; 

ASK  METHOD  GiveMeanRCS(IN  ThisThreat :  STRING)  REAL 
END  OBJECT; 

EnvironmentObj  =  OBJECT 

RandTime  :  REAL; 

ClutterVar  ;  REAL; 

SystemLoss, 

AtmosLoss  :  REAL; 

ActiveThreat, 

ActivePlat  ;  PlatformObj; 

FutureLocale  :  LocationRecType; 

MyRand  :  RandomObj; 

ASK  METHOD  Objinit; 

ASK  METHOD  GiveRCS(IN  NewThreat ;  STRING);  REAL 
ASK  METHOD  Clutter()  ;  REAL; 

ASK  METHOD  SysLossQ  :  REAL; 

ASK  METHOD  AtmLoss()  :  REAL; 

ASK  METHOD  SetCurrentPlats(IN  P  ;  PlatformObj;  IN  T  ;  PlatformObj); 

ASK  METHOD  SetEntryLocation(IN  Locale  :  LocationRecType); 

ASK  METHOD  FindDetectTime(IN  Pulseint :  INTEGER;IN  ThiaPlat ;  STRING; 

IN  ThisSensor :  STRING; 

IN  ThisTarget :  PlatformObj;  IN  Ro  :  REAL); 
ASK  METHOD  PassDetectionRCS(IN  ThisThreat .  PlatformObj)  REAL- 
END  OBJECT; 

RadarObj  =  OBJECT 

Yb  ;REAL; 

ReceiverTemp  :  REAL; 

Rmax  :  ARRAY  INTEGER  OF  REAL;  {MAX  RANGE  FOR  EACH  PLATFORM  TYPE} 
SNRatio  ;  ARRAY  INTEGER  OF  REAL; 

RCSStatKeeper ;  StatObj; 

RandVarl  :  RandomObj; 

ASK  METHOD  Objinit; 

ASK  METHOD  Gain(IN  ThisSensor  :  SensorlnfoRecType)  :  REAL; 

ASK  METHOD  FindRo(IN  TgtType  ;  STRING;IN  SensorTypeName  :  STRING)  :  REAL; 
ASK  METHOD  TfaermalVar(IN  Csensor  :  SensorlnfoRecType)  ;  REAL; 

ASK  METHOD  Interference(IN  CSensor ;  SensorlnfoRecType) :  REAL; 

ASK  METHOD  FmdBiasLevel(IN  b  ;  INTEGER;  IN  NPrime  ;  REAL)  :  REAL; 

ASK  METHOD  Integrate Yb(IN  fofYb  :  REAL;  IN  N  ;  INTEGER)  :  REAL; 

ASK  METHOD  ForElevationAngle{)  ;  REAL; 

ASK  METHOD  GetNoise(IN  noise  ;  REAL)  :  REAL; 

ASK  METHOD  FmdTheSNRatio(IN  TP  :  LocationRecType;  IN  PP  ;  LocationRecType; 

IN  SensParam :  SensorlnfoRecType; 

IN  TgTObj  :  PlatformObj); 


END  OBJECT; 

VAR 

Environment :  EnvironmentObj; 


{  Need  to  be  newed  some;vhere} 
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Radar  ;  RadarObj;  {  Need  to  be  newed  somewhere} 

Threat  .  ThreatObj;  { Need  to  be  newed  somewhere} 

Threshold  ;  REAL; 
k  :REAL; 

END  MODULE. 


Implementation  module  lAARONl.mod: 


IMPLEMENTATION  MODULE  AARONl; 

FROM  GM  IMPORT  LocationRecType; 

FROM  Plat  IMPORT  PlatformObj; 

FROM  MathMod  IMPORT  LOGl  0,  POWER,  EXP,  LN,  pi,  COS,  SQRT; 

{  use  this  to  get  info  from  the  platform  records} 

FROM  PList  IMPORT  MasterPlatformTnfnT  ict  Pl!>ri7-vr»vi¥»iiv.T  wrtu; 


FROM  PList  IMPORT  MasterPlatformlnfoList,  PlatfonnInfoListObj, 

PlatformInfoRecType; 

FROM  DetRng  IMPORT  AttributesRecType,  MasterSensorList,  SensorMasterListObj, 

SensorInfoRecType,  DetectionRangePairObj; 

FROM  MathFunc  IMPORT  MyMath; 

FROM  Write  IMPORT  WriteLine,  WriteData; 

FROMAARONCDF  IMPORT  BlackBox; 

FROM  AARONTIME  IMPORT  MasterEntryTime; 

FROM  TypeLis  IMPORT  MasterPlatformTypeList,  PlatformTypeRecType- 

FROM  Space  IMPORT  Sensor; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  Manuv  IMPORT  NamedManeuverObj,  ManeuverListObj,  NamedPathObj, 

ManeuveringPlatformObj; 

FROM  SimpleStats  IMPORT  StatObj; 

FROM  SimExp  IMPORT  ExperimentManager; 


FROM  Space 
FROM  RandMod 
FROM  Manuv 

FROM  SimpleStats 
FROM  SimExp 


OBJECT  ThreatObj; 
ASK  METHOD  Objinit; 
BEGIN 

END  METHOD; 


ASK  METHOD  GiveMeanRCS(lN  ThisThreat :  STRING)  REAL 
VAR 

Rec ;  PlatformInfoRecType; 

BEGIN 

Rec  ;=  ASK  MasterPlatformlnfoList  TO  GivePlatform(ThisThreat)- 
RETURN(Rec.RCS); 

END  METHOD; 

END  OBJECT; 


OBJECT  EnvironmentObj; 

ASK  METHOD  Objinit; 

BEGIN 

ClutterVar  :=  0.0; 

SystemLoss  ;=  1.26; {Says  system  losses,  IdB} 
AtmosLoss  :=  3.16;{Says  atmos  losses,  5dB} 
NEW(MyRand);  (Used  for  the  RCS} 
RandTime  :=  0.0;  {Initialized  to  zero} 

END  METHOD; 
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ASK  METHOD  GiveRCS(IN  NewThreat :  STRING);  REAL- 
BEGIN 

RETURN(ASK  Threat  GiveMeanRCS(NewThreat)) 

END  METHOD; 

ASK  METHOD  ClutterQ  :  REAL; 

BEGIN 

RETURN(ClutterVar); 

END  METHOD; 

ASK  METHOD  SysLossO  :  REAL; 

BEGIN 

RETURN(SysteniLoss); 

END  METHOD; 

ASK  METHOD  AtmLossO  :  REAL; 

BEGIN 

RETURN(AtmosLoss); 

END  METHOD; 

ASK  METHOD  SetCurrentPIats(IN  P  :  PlatformObj;  IN  T  ;  PlatformObj); 
{This  method  sets  the  active  participants  platform  objects) 

BEGIN 

ActiveThreat  :=  T; 

ActivePlat  :=  P; 

END  METHOD; 

ASK  METHOD  SetEntryLocation(lN  Locale  ;  LocationRecType)' 

BEGIN 

FutureLocale  :=  Locale; 

END  METHOD; 


ASK  METHOD  FindDetectTime(IN  Pulseint :  INTEGER;IN  ThisPlatID  :  STRING; 

IN  ThisSensor  :  STRING;IN  ThisTarget :  PlatformObj; 
IN  Ro:  REAL); 

VAR 


Time,  SlantRangeTime, 

speed 

temprec 

IDRec 

TrueRange, 

RadianAngle 

Xnot,  Znot,  a,b,c,VSI 

MyPosit, 

ThreatPosifThreatVel 

VelocityZ,VX 


;REAL; 

:  PlatformTypeRecType; 
:  PlatformlirfoRecType; 

REAL; 

:REAL; 

:  LocationRecType; 
:REAL; 


ThisPath 

TheseManeuvers 

NumberOMan 

SensorParameters 

BEGIN 


;  NamedPathObj; 

:  ManeuverListObj; 

:  INTEGER; 

:  SensorInfoRecType; 


SensorParameters  :=  ASK  MasterSensorList  GiveSensor(ThisSensor); 

ASK  BlackBox  TO  SetLookTimes(SensorParameters, {TheseManeuvers,}  Ro); 

SlantRangeTime  :=  ASK  BlackBox  TO  FindARange(PulseInt);  {Returns  random  range  of  detection) 
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IDRec  ;=  ASK  MasterPlatformlnfoList  GivePlatform(ThisPlatID); 

{Pass  Time  to  EDetect.mod} 

RandTime  :=  SlantRangeTime; 

{Store  times  here} 

ASK  MasterEntiyTime  TO  AddRecord(IDRec,PIatformName,ThisSensor, 

ASK  ThisTarget  ID,  SlantRangeTime); 


END  METHOD; 

ASK  METHOD  PassDetectionRCS(IN  ThisThreat :  PlatformObj) :  REAL" 

VAR 

muRCS,  R  :  REAL; 

RCSatDetection :  REAL; 

BEGIN 

muRCS  :=  ASK  SELF  GiveRCS(ASK  ThisThreat  ED); 

R  :=  ASK  MyRand  UniformReal(0.0,l .0); 

RCSatDetection  :=  -muRCS  *  LN(R); 

RETURN(RCSatDetection); 

END  METHOD; 

END  OBJECT; 

OBJECT  RadaiObj; 

ASK  METHOD  Objinit; 

BEGIN 

k  :=  1.38/POWER(10.0,23.0); 

NE  W(RC  S  StatKeeper) ; 

NEW(RandVarl); 

ASK  RandVarl  TO  SetSeed(21 16429); 

END  METHOD; 

ASK  METHOD  ThermalVar(IN  Csensor  :  SensorMoRecType)  :  REAL; 

VAR 

value :  REAL; 

BEGIN 

value  :=  k*Csensor.Attributes.ReceiverTemp*Csensor. Attributes  FilterBW  • 
RETURN(value); 

END  METHOD; 

ASK  METHOD  Interference(IN  CSensor  ;  SensorMoRecType)  •  REAL 
VAR 

IntVar :  REAL; 

BEGIN 

IntVar  :=  ASK  SELF  Thermal Var(CSensor)  +  ASK  Environment  TO  ClutterQ' 
RETURN(IntVar); 

END  METHOD; 

ASK  METHOD  FindRo(IN  TgtType :  STRING; 

IN  SensorTypeName  :  STRING):  REAL; 

VAR 

sigma,  terml,  term2,  term3,A :  REAL; 
i  :  INTEGER; 

CurrentSensor  :  SensorMoRecType; 

Noise,  Clutter  ;  REAL; 

B  :  INTEGER; 
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BEGIN 

CurrentSensor  ;=  ASK  MasterSensorList  TO  GiveSensor(SensorT5^eName); 
sigma  ;=  ASK  Environment  TO  GiveRCS(TgtType); 

terml  ;=  CurrentSensor.Attributes.AveragePower  *  CurrentSensor. Attributes.Aperture  * 

ASK  SELF  Gain(CurrentSensor)*sigma  *  CurrentSensor.Attributes.EfBciency  *  0.5  * 
FLOAT(CurrentSensor.Attributes.PulseInt); 
term2  :=  (16.0*POWER(pi,FLOAT(2)))  *  Interference(CurrentSensor) ; 
term3  :=  ASK  Environment  SysLoss()  *  ASK  Environment  AtmLoss(); 

B  :=  CurrentSensor. Attributes.PulseInt; 

A  :=  CurrentSensor.  Attributes.FalseAlarmNumber; 

Yb  ;=  ASK  SELF  FindBiasLevel(B,  A);  (Sets  the  Threshold} 

(Returns  the  Max  det  range  based  on  RCS  and  radar  parameters} 
RETURN(POWER(terml/(term2  *  term3),0.25)  / 1852.0);  (Convert  from  meters} 

END  METHOD; 

ASK  METHOD  Gain(IN  ThisSensor  :  SensorInfoRecType):  REAL; 

(Need  to  pass  the  actual  record} 

VAR 

ap  :  REAL; 
lambda ;  REAL; 

BEGIN 

ap  :=  ThisSensor. Attributes.Aperture; 
lambda  :=  ThisSensor. Attributes.WaveLength; 

RETURN{(4.0  *  pi  *  ap)/POWER{lambda,FLOAT(2))); 

END  METHOD; 

ASK  METHOD  FindBiasLevel(IN  b  ;  INTEGER;  IN  NPrime  ;  REAL)  :  REAL; 

(  NPrime=  FALSE  ALARM  NUMBER,  b=  PULSES  TO  INTEGRATE} 

VAR 

numerator,  denoml,  denom2  :  REAL; 
a,  C  ;  REAL; 

BEGIN 

a  ;=  POWER(0.5,  (1 .0/NPiime));  (This  is  the  value  to  check  against} 

(a  is  FofYb} 

C  ;=  ASK  SELF  TO  Integrate Yb{a,  b);(Pass  FofYb  and  Pulses  to  Integrate} 

RETURN(C); 

END  METHOD; 

ASK  METHOD  IntegrateYb(IN  fofYb  ;  REAL;  IN  N  :  INTEGER)  :  REAL; 

(Integrate  For  the  purpose  of  finding  the  Threshold} 


VAR 

numerator,  denoml,  denom2  :  REAL; 

a,  b,  h,  temp,  sum,  X,  dumy  ;  REAL; 

answemew,  y,  delta  :  REAL; 

FoimdSolution  :  BOOLEAN; 

littlen,  i  :INTEGER; 

BEGIN 

littlen  :=  800; 


FoundSoIution  ;=  FALSE; 
a:=  0.0; 
b:=  0.1; 

h:=(b-a)/FLOAT(littlen); 

temp  ;=  (EXP(-1.0*b))  *  (POWER(b,  FLOAT(N  -  1)));  (sum  of  f(a)  and  f(b)} 
sum  ;=  ASK  MyMath  TO  ComputeFactorial(N  - 1,  temp);  (this  is  the  last  term} 
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FOR  i;=  1  TO  (littlen  -  1) 

X:=h*FLOAT(i); 

dumy:=  (EXP(-1.0  *  X))  *  {POWER(X,  FLOAT(N-l)))  ;{2  times  y} 
y ASK  MyMath  TO  ComputeFactoria](N  -  l,dumy); 
sum:=  sum  +  (2.0  *  y); 

END  FOR; 

answemew  :=  (h/2.0)  *  sum; 

delta  ;=  0.1;  {  This  is  the  increase  in  Yb  until  the  right  value  is  found  1 
REPEAT 

b  :=  b  +  delta;  . 
h:=(b-a)/FLOAT(littlen); 

temp  ;=  (EXP(-1.0*b))  *  (POWER(b,FLOAT(N-l)));  (sum  of  f(aj  and  f(b)} 
sum  ASK  MyMath  TO  ComputeFactorial(N  -  1  temp)- 
FORi;=  1  TO{Uttlen- 1) 

X:=h*FLOAT(i); 

dumy;=  (EXP(-1.0*X))  *  (POWER(X,FLOAT(N-l)))  ;{2  times  y) 
y  :=  ASK  MyMath  TO  ComputeFactorial(N  -  l,dumy); 
sum;=  sum  +  (2.0  *  y); 

END  FOR; 

answemew  ;=  (h/2.0)  *  sum; 

IF  fofVfb  <  answemew 
Threshold  :=  b; 

FoundSolution  :=  TRUE; 

RETURN(b);  {This  b  is  the  thresholded  value} 

END  IF; 

UNTIL  FoundSolution; 

END  METHOD; 

ASK  METHOD  ForElevationAngleQ  :  REAL 
VAR 

Angle  .  REAL; 

BEGIN 

NEW(Sensor); 

Angle  ;=  ASK  Sensor  TO  GiveElevationAngle() 

DISPOSE(Sensor); 

RETURN(Angle); 

El^  METHOD; 

ASK  METHOD  GetNoise(IN  noise  :  REAL)  :  REAL; 

{Uses  the  invers  transform  of  the  Weibull} 

VAR 
U ,  Beta, 

RandWeibull :  REAL; 

BEGIN 

U  :=  ASK  RandVarl  UniformReal(0.0, 1 .0); 

Beta  :=  SQRT(2.0*noise); 

RandWeibull  :=  Beta*POWER(-LN(U),0.5); 

RETURN(RandWeibull); 

END  METHOD; 

ASK  METHOD  FindTheSNRatio(IN  TP  :  LocationRecType;  IN  PP  ;  LocationRecType; 

IN  SensParam  :  SensorInfoRecType; 

IN  TgTObj  :  PlatformObj); 
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{For  analysis  purposes  find  the  S/N  Ratio) 

VAR 

term  I,  term2, 
terms,  term4, 

SNRatio  :  REAL; 

TheRange  :  REAL; 

DetRCS  :  REAL; 

TheNoise, 

TheSignaI,SplusN  :  REAL; 

BEGIN 

TheRange  ;=  SQRT(POWER(TP.x-PP.x,2.0)  +  POWER(TP.y-PP.y,2.0)  + 

POWER(TP.z-PP.z,2.0))  *  1852.0;  {Convert  units  to  meter) 

DetRCS  :=  ASK  Environment  TO  PassDetectionRCS(TgTObj);{Here  use  rand  RCS) 
terml  ;=  SensParam.Attributes.AveragePower*Gain(SensParam)*DetRCS* 
SensParam.Attributes.Aperture*SensParam.Attributes.EfBciency* 
FLOAT(SensParam.Attributes.PulseInt); 
term2  :=POWER((4.0*pi*POWER(TheRange,2.0)),2.0); 
terms  :=  Interference(SensParam); 

term4  ;=  ASK  Environment  SysLossQ  *  ASK  Enviromnent  AtmLossO; 

SNRatio  ;=  terml/(term2  *  termS  *  tenn4); 

TheNoise  :=  GetNoise(termS); 

TheSignal  :=tennl/(term2  ♦  term4); 

SplusN  :=  (TheSignal  +  TheNoise)*POWER(10.0,14.0); 

WriteData(TheSignal*PO\VER{l  0.0, 14.0)); 

ASKRCSStatKeeperTO  GetSampl^TheSignal*POWER(10.0,14.0));  {Grab  Stats) 

IF  ((ASK  RCSStatKeeper  HWQ  /  ASK  RCSStatKeeper  Mean())  <= 

ASK  ExperimentManager  DesiredPrecision)  AND  (ASK  RCSStatKeeper  N  >  1) 
ASK  RCSStatKeeper  TO  Output; 

ASK  ExperimentManager  TO  SetBoolean; 

END  IF; 

END  METHOD; 

END  OBJECT; 

END  MODULE. 

lAARONCDF.mod  builds  the  cumulative  distribution  function  from  which  the 
inverse  transform  method  determines  the  random  time  of  detection. 

Definition  module  DAARONCDF.mod: 

DEFINITION  MODULE  AARONCDF; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  Manuv  IMPORT  ManeuverListObj; 

FROM  DetRng  IMPORT  AttributesRecType,SensorInfoRecType; 

FROM  GM  IMPORT  LocationRecType; 

{  THE  PURPOSE  OF  THIS  MODULE  IS  TO  GENERATE  THE  RV  Time) 

TYPE 

CDFRec  =  RECORD 

Range  :  REAL; 

NoDetect :  REAL; 

CumPd  :  REAL; 
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PDF  :  REAL; 

END  RECORD; 

CDFObj  =  OBJECT 

Pd  :  ARRAY  INTEGER  OF  REAL; 

Ufour  :  ARRAY  INIEGER  OF  REAL; 

CumProb  :  ARRAY  INTEGER  OF  CDFRec; 

Ro  :  REAL; 

RandVar  :  RandomObj; 

Ratio  :  ARRAY  INTEGER  OF  REAL; 

LookRanges  :  ARRAY  INTEGER  OF  REAL; 

TimeToRange  ;  ARRAY  INTEGER  OF  REAL; 

K  :  INTEGER; 

ASK  METHOD  Objinit; 

ASK  METHOD  FmdPosTime(IN  CSpeed  :  LocationRecType; 

IN  TPos  ;  LocationRecType; 

IN  PPos  :  LocationRecType; 
INRng  ;  REAL);  REAL; 

ASK  METHOD  SetLookTimes(IN  sensorrec  :  SensorInfoRecType; 

IN  NewRo  ;  REAL  ); 

ASK  METHOD  FindARange(IN  N  :  INTEGER)  ;  REAL; 

ASK  METHOD  FmdSingleHitPd(IN  N :  INTEGER)  :  REAL; 

ASK  METHOD  BuildCDF(IN  N  ;  INTEGER)  :  REAL; 

ASK  METHOD  FindRandRange :  REAL; 

END  OBJECT; 

CONST  Partitions  =  900;  {Set  intervals  for  max  number  of  looks) 

VAR 

BlackBox :  CDFObj; 

END  MODULE. 

Implementation  module  lAARONCDF.mod: 

IMPLEMENTATION  MODULE  AARONCDF; 

FROM  MathMod  IMPORT  POWER,  FLOOR,  EXP,  SQRT,  SIN,  COS,  pi,  e; 

FROM  AARON  1  IMPORT  Environment,  Threshold; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  Manuv  IMPORT  ManeuverListObj,  NamedManeuverObj; 

FROM  DetRng  IMPORT  AttributesRecType,SensorInfoRecType; 

FROM  GM  IMPORT  LocationRecT)fpe; 

FROM  AngleU  IMPORT  RadToDeg,  RadEastToDegNorth,  BearingToLocation; 
FROMMathFunc  IMPORT  MyMath; 

FROM  PList  IMPORT  MasterPIatformlnfoList,  PlatformlnfoListObj, 

PlatformlnfoRecType; 

OBJECT  CDFObj; 

ASK  METHOD  Objinit; 

VAR 

i:  INTEGER; 

BEGIN 

NEW(RandVar); 

ASK  RandVar  TO  SetSeed(21 16429302); 

END  METHOD; 
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ASK  METHOD  FindPosTime(  IN  RSpeed  :  LocationRecType;  {Relative  Speed} 

IN  TPos  ;  LocationRecType; 

IN  PPos  :  LocationRecType; 

INRng  :  REAL);  REAL; 

VAR 

a,  b,  c,  time  :  REAL; 

BEGIN 

a  ;=  POWER(RSpeed.x,2.0)  +  POWER(RSpeed.y,2.0)  +  POWER(RSpeed.z,2.0); 
b  :=  2.0*((TPos.x*RSpeed.x)  +  (TPos.y*RSpeed.y)  +  (TPos.z*RSpeed.z)  + 

{-PPos.x*RSpeed.x)  +  (-PPos.y*RSpeed.y)  +  (-PPos.z*RSpeed.z)); 
c  :=  POWER(TPos.x,2.0)  +  POWER(TPos.y,2.0)  +  POWER(TPos.z,2.0)  + 
PO\VER(PPos.x,2.0)  +  POWER^Pos.y,2.0)  +  POWER^Pos.z,2.0)  - 
2.0*((TPos.x*PPos.x)  +  (TPos.y*PPos.y)  +  (TPos.z*PPos.z))  -  POWER(Rng,2.0); 
time  ;=  ASK  MyMath  Quad(a,b,c); 

RETlJRN(tune); 

END  METHOD; 

ASK  METHOD  SetLookTimes(IN  sensorrec ;  SensorlnfoRecType; 

IN  NewRo  ;  REAL  ); 

VAR 

i,j  INTEGER; 


curRo,  nextRange, 

ElaspedTime,  TgtCourse, 

PCourse,  RangeCheck, 

X,  y,  2,  t,  V 

velocity,  curBigXY, 
curBigZ,  DestBigXY, 
DestBigZ,  Bearing, 
NextBearing,  speedZ, 

:REAL; 

speedXY,  Timeint 
curmanuv, 

;REAL; 

nextmanuv,  pnmanuv 

PlatCenter, 

curTgtPosit, 

:  NamedManeuverObj; 

DestTgtPosit,  TDestTgtPosit, 

EPTgt,  IPPlat, 

MaxTgtPosit, 

PlatSpeed,  future, 

SpeedVec,  TgtSpeed, 

RoPosit ,  PDest 
Properlnterval, 

:  LocationRecType; 

Inbound 

:  BOOLEAN; 

counter 

:  INTEGER; 

TheseAttributes 

:  AttributesRecType; 

TgtManuv,  Pmanuv 

:  ManeuverListObj; 

Ttemprec, Ptemprec 
BEGIN 

:  PlatformlnfoRecType; 

TgtManuv  :=  ASK  Environment  ActiveThreat.CurrentPath  ManeuverList; 

Ttemprec  :=  ASK  MasterPiatformInfoList  GivePlatform(ASK  Environment.  ActiveThreat  ID)- 
NEW(IPTgt); 

IPTgt.x  :=  Ttemprec.InitiaILocation.x; 

IPTgt.y  ;=  Ttemprec.InitialLocation  y; 

IPTgt.z  :=  Ttemprec.InitiaILocation.z; 
counter  :=  0; 
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Ro  ;=  NewRo; 

Timeint  :=  sensorrec.Attributes.ScanTime;{Time  allowed  for  movement  between  scans} 
curmanuv  :=  ASK  TgtManuv  FirstO;{  Always  at  least  1  maneuver} 
curTgtPosit  ;=  IPTgt; 

nextmanuv  :=  ASK  TgtManuv  Next(cunnanuv); 

{  ***  Plat  is  not  moving  so  find  its  first  speed} 

NEW(PlatSpeed); 

Ptemprec  :=  ASK  MasterPlatformlnfoList  GivePlatfonn(ASK  Environment. ActivePlat  ID); 
PlatCenter  ;=  Ptemprec.InitialLocation; 

Pmanuv  :=  ASK  Environment.  ActivePlat.  CurrentPath  ManeuverList; 
pnmanuv  :=  ASK  Pmanuv  FirstO; 

PDest  :=  ASK  pnmanuv  DestinationForStraightShot; 

PCourse  :=  BearingToLocation(PlatCenter,  PDest); 

PlatSpeed.z  ;=  (ASK  pnmanuv  ClimbDiveRate)/6000.0; 
velocity  :=  ASK  pnmanuv  NewSpeed; 

speedXY  ~  SQRT(POWER(velocity,2.0)  -  POWER(speedZ,2.0)); 

PlatSpeed.x  ;=  speedXY  *  SIN(PCourse  *  pi/180.6); 

PlatSpeed.y  :=  speedXY  *  COS(PCourse  *  pi/180.0); 

Inbound  :=  TRUE; 

Bearing  :=  BearingToLocation(PlatCenter,  IPTgt); {Assumes  Straight  Shot} 

Properinterval  ;=  FALSE; 

NEW(future); 

NEW(SpeedVec); 

NEW(TgtSpeed); 

NEW(RoPosit); 

NEW(DestTgtPosit); 

NEW^ookRanges,  1 .  .Partitions); 

NEW(TimeToRange,  1 .  .Partitions); 

{  ***  Find  the  correct  starting  interval} 

WHILE  (curmanuv  o  NILOBJ)  AND  (Inbound) 

{ Sets  the  climb/descent  speed  and  converts  to  Mile  per  min,  ML  in  man.dat  doesn't  work} 
curBigXY  :=  SQRT(POWER((curTgtPosit.x  -  PlatCenter.x),2.0)  + 

POWER((curTgtPosit.y  -  PlatCenter. y),2.0)); 
curBigZ  ;=  curTgtPosit.z  -  PlatCenter.z ; 

{  ****  This  is  for  true  destination  } 

TDestTgtPosit  :=  ASK  curmanuv  DestinationForStraightShot; 

{  ***Sets  the  inbound  speed  Vector  of  the  target} 

TgtCourse  :=  BearingToLocation(curTgtPosit,  TDestTgtPosit); 

TgtSpeed.z  :=  (ASK  curmanuv  ClimbDiveRate)/6000.0; 
velocity  ;=  ASK  curmanuv  NewSpeed; 

speedXY  :=  SQRT(POWER(velocity,2.0)  -  POWER(speedZ,2.0)); 

TgtSpeed.x  :=  speedXY  *  SIN(TgtCourse  *  pi/180.0); 

TgtSpeed.y  :=  speedXY  *  COS(TgtCourse  *  pi/180.0); 

{  ***  Sets  the  Rate  of  Closure  Vector } 

SpeedVec.x  :=  TgtSpeed.x  -  PlatSpeed.x; 

SpeedVec.y  :=  TgtSpeed.y  -  PlatSpeed.y; 

SpeedVec.z  :=  TgtSpeed.z  -  PlatSpeed.z; 

t 

{  Rework  and  solve  for  relative  destination  and  check,  v  is  the  time  to  get  to 

true  target  destination  from  current  TgtPosit} 
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V  :=  SQRT(POWER(TDestTgtPosit.x  -  curTgtPosit.x,2.0)  + 

POWER(TDestTgtPosit.y  -  curTgtPosit.y,2.0)  + 

POWER(TDestTgtPosit.y  -  curTgtPosh  y,2.0))  /  velocity; 

IF  nextmanuv  o  NILOBJ 

DestTgtPosit.x  :=  curTgtPosit.x  +  (SpeedVec.x*v); 

DestTgtPosit.y  :=  curTgtPosit.y  +  (SpeedVec.y*v); 

DestTgtPosit.z  :=  curTgtPosit.z  +  (SpeedVec.z*v); 

DestBigXY  :=  SQRT(POWER((DestTgtPosit.x  -  PlatCenter.x),2.0)  + 

POWER(^estTgtPosit.y  -  PlatCenter.y),2.0)); 

DestBigZ  :=  DestTgtPosit.z  -  PlatCenter.z; 

ELSE 

DestTgtPosit.x  :=  TDestTgtPosit.x; 

DestTgtPosit.y  ;=  TDestTgtPosit.y ; 

DestTgtPosit.z  :=  TDestTgtPosit.z; 

DestBigXY  :=  SQRT(POWER{(DestTgtPosit.x  -  PlatCenter.x),2.0)  + 

POWER(^estTgtPosit.y  -  PlatCenter.y),2.0)), 

DestBigZ  :=  DestTgtPosit.z  -  PlatCenter.z; 

END  IF; 

IF  {SQRT(POWER(curBigXY,2.0)  +  POWER(curBigZ,2.0))  <= 

SQRT(POWER(DestBigXY,2.0)  +  POWER(DestBigZ,2.0)))  AND  (Inbound) 
MaxTgtPosit  ;=  curTgtPosit; 

{  ***  t  is  the  time  to  get  to  Ro  from  MaxTgtPosit) 
t  :=  FindPosTime(SpeedVec,  MaxTgtPosit,  PlatCenter,  Ro); 

{  ***  Find  the  coordinates  of  Ro  } 

RoPoat.x  :=  MaxTgtPosit.x  +  (SpeedVec.x*t); 

RoPosit:y  ;=  MaxTgtPosit.y  +  (SpeedVec.y*t); 

RoPosit.z  :=  MaxTgtPosit.z  +  (SpeedVec.z*t); 

RangeCheck  :=  SQRT(POWER(curBigXY,2.0)  +  POWER(curBigZ,2.0)); 

Properinterval  :=  TRUE; 

REPEAT 

future.x  ;=  curTgtPosit.x  +  (Sj>eedVec.x*TimeInt);(New  x  posit) 
future.y  :=  curTgtPosit.y  +  (SpeedVec.y*TimeInt);{New  y  posit) 
foture.z  ;=  curTgtPositz  +  (SpeedVec.z*TimeInt);{New  z  posit) 
nextRange  ;=  SQRT(POWER(fiiture.x  -  PlatCenter.x,2.0)  + 

POWER(future.y  -  PlatCenter.  y,2.0)  + 

POWER(future.z  -  PlatCenter.z,2.0)); 
curTgtPosit  ;=  future; 

{  Check  if  the  range  is  decreasing,  if  so,  add  to  array  if  R  <  Ro) 

IF  nextRange  <=  RangeCheck 

curBigXY  :=  SQRT(POWER((curTgtPosit.x  -  PiatCenter.x),2.0)  + 
POWER((curTgtPosit.y  -  PlatCenter. y),2.0)); 
curBigZ  :=  curTgtPosit.z  -  PlatCenter.z; 

RangeCheck  :=  nextRange; 

IF  nextRange  <=  Ro 
counter  :=  counter  +  1; 

LookRanges[counter]  :=  nextRange; 

{  ***  Time  to  go  from  Ro  to  all  other  accepted  positions  } 

IF  counter  o  1 

TimeToRange[counter]  :=  (SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  PO\VER(future.y  -  RoPosit.y,2.0)  + 
POWER(future.z  -  RoPosit.z,2.0)) 

/  (SQRT(POWER(SpeedVec.x  ,2.0)  + 
POWER(SpeedVec.y,2.0)  + 
POWER(SpeedVec.z,2.0))))  + 
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TimeToRange[counter  -  1] ; 

RoPosit.x  :=  future.x; 

RoPosit.y  :=  fiiture.y; 

RoPosit.z  :=  future.z; 

ELSE 

{  ***  Time  to  go  from  Ro  to  first  accepted  position} 

TimeToRange[counter]  ;=  SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  POWER(future.y  -  RoPosit.y,2.0) 

+  POWER^future.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 
POWER(SpeedVec.y,2.0)  + 
POWER{SpeedVec,z,2.0)); 

RoPosit.x  :=  future.x; 

RoPosit.y  ;=  future.y; 

RoPosit.z  :=  future.z; 

END  IF; 

END  IF; 

ELSE 

Inbound  ;=  FALSE;{Tgt  passed  ovhd) 
cunnanuv  ;=  nextmanuv; 

END  IF; 

UNTIL  NOT(Inbound); 

(Force  the  target  outbound  after  min  range} 
inbound  :=  FALSE;{Tgt  passed  ovhd} 

{  Removed  while  loop  because  since  in  this  manuever  the  target  will  pass 

ovhd  you  do  not  want  to  consider  a  next  manuever} 
(for  when  Max  is  less  than  Min,  the  normal  configuration  execute  the  following} 

ELSIF  (Ro  >  SQRT(POWER(DestBigXY,2.0)  +  POWER(DestBigZ,2.0)))  AND  (Inbound) 
MaxTgtPosit  ;=  curTgtPosit; 

{  ***  t  is  the  time  to  get  to  Ro  from  MaxTgtPosit} 
t  :=  FmdPosTime(SpeedVec,  MaxTgtPosit,  PlatCenter,  Ro); 

{  ***  Find  the  coordinates  of  Ro } 

RoPosit.x  :=  MaxTgtPosit.x  +  (SpeedVec.x*t); 

RoPosit.y  :=  MaxTgtPosit.y  +  (SpeedVec.y*t); 

RoPosit.z  :=  MaxTgtPosit. z  +  (SpeedVec.z*t); 

RangeCheck  :=  SQRT(POWER(curBigXY,2.0)  +  POWER(curBigZ,2.0)); 

Properinterval  :=  TRUE; 

REPEAT 

fiiture.x  :=  curTgtPosit. x  +  (SpeedVec.x*TimeInt);{New  x  posit} 
future.y  :=  curTgtPosit.y  +  (SpeedVec.y*TimeInt);{New  y  posit} 
fijture.z  :=  curTgtPosit.z  +  (SpeedVec.z*TimeInt);{New  z  poat} 
nextRange  :=  SQRT(POWER(future.x  -  PlatCenter.x,2.0)  + 

POWER(future.y  -  PlatCenter.y,2.0)  + 

POWER(future.z  -  PlatCenter.z,2.0)); 
curTgtPosit  :=  future; 

{  ***  Check  if  the  range  is  decreasing,  if  so,  add  to  array  if  R  <  Ro} 

IF  nextRange  <=  RangeCheck 

curBigXY  :=  SQRT(POWER((curTgtPosit.x  -  PlatCenter .x),2.0)  + 

POWER((curTgtPosit.y  -  PlatCenter. y),2.0)); 
curBigZ  ;=  curTgtPosit.z  -  PlatCenter.z; 

RangeCheck  ;=  nextRange; 
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IF  nextRange  <=  Ro 
counter  ;=  counter  +  I; 

LookRanges[counter]  ~  nextRange; 

{  ***  Time  to  go  from  Ro  to  all  other  accepted  positions  } 

IF  counter  o  1 

TimeToRangefcounter]  :=  (SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  PO\VER(future.y  -  RoPosit.y,2.0) 

+  PO\VER{future.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 

POWER(SpeedVec.y,2.0)  + 
POWER(SpeedVec.z,2.0))) 

+  TimeToRangefcounter  -  1] ; 

RoPosit.x  ;=  fiiture.x; 

RoPosit.y  :=  fiiture.y; 

RoPosit.z  :=  fiiture.z; 

ELSE 

{  ***  Time  to  go  from  Ro  to  first  accepted  position} 

TimeToRangefcounter]  :=  SQRT(PO\\^R(fiiture.x  -  RoPosit.x,2.0) 

+  POWER(future.y  -  RoPosit.y,2.0) 

+  POWER(future.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 

POWER(SpeedVec.y,2.0)  + 

POWER(SpeedVec.z,2.0)); 

RoPosit.x  :=  fiiture.x; 

RoPosit.y  :=  fiiture.y; 

RoPosit.z  :=  fiiture.z; 

END  IF; 

END  IF; 

ELSE 

Inbound  ;=  FALSE;  (Tgt  passed  ovhd} 
curmanuv  :=  nextmanuv; 

END  IF; 

UNTIL  (nextRange  <  SQRT(POWER(DestBigXY,2.0)  +  POWER(DestBigZ,2.0)))  OR  (NOT(Inbound)); 
WHILE  nextmanuv  o  NDLOBJ  {Move  to  next  Maneuver} 

MaxTgtPosit curTgtPosit; 

curmanuv  :=  nextmanuv; 

nextmanuv  ;=  ASK  TgtManuv  Next(curmanuv); 

{  ****  This  is  for  true  destination  } 

TDestTgtPosit  ;=  ASK  curmanuv  DestinationForStraightShot; 

TgtCourse  :=  BeaiingToLocation(curTgtPosit,  TDestTgtPosit); 

TgtSpeed.z  :=  (ASK  curmanuv  ClimbDiveRate)/6000.0; 
velocity  ;=  ASK  curmanuv  NewSpeed; 

speedXY  :=  SQRT(POWER(velocity,2.0)  -  POWER(speedZ,2.0)); 

TgtSpeed.x  ;=  speedXY  *  SIN(TgtCourse  *  pi/180.0); 

TgtSpeed.y  :=  speedXY  *  COS(TgtCourse  *  pi/180.0); 

SpeedVec.x  :=  TgtSpeed.x  -  PlatSpeed.x; 

SpeedVec.y  :=  TgtSpeed.y  -  PlatSpeed.y; 

SpeedVec.z  :=  TgtSpeed.z  -  PlatSpeed.z; 

V  :=  SQRT(POWER(TDestTgtPosit.x  -  curTgtPosit.x,2.0)  + 

POWER(TDestTgtPosit.y  -  curTgtPosit. y,2.0)  + 

POWER(TDestTgtPosit.y  -  curTgtPosit.y,2.0))  /  velocity; 

IF  nextmanuv  o  NILOBJ  > 

DestTgtPosit.x  :=  curTgtPosit.x  +  (SpeedVec.x*v); 
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DestTgtPosit.y  ;=  curTgtPosit  y  +  (SpeedVec.y*v); 

DestTgtPosit.z  :=  curTgtPosit.z  +  (SpeedVec.z*v); 

DestBigXY  :=  SQRT{POWER((DestTgtPosit.x  -  PIatCenter.x),2.0)  + 
POWER((PestTgtPosit.y  -  PlatCenter.y),2.0)); 

DestBigZ  :=  DestTgtPosit.z  -  PlatCenter.z; 

ELSE 

DestTgtPosit.x  ;=  TDestTgtPosit.x; 

DestTgtPosit.y  ;=  TDestTgtPosit.y ; 

DestTgtPosit.z  :=  TDestTgtPosit.z; 

DestBigXY  :=  SQRT(POWER((DestTgtPosit.x  -  PlatCenter.x),2.0)  + 
POWER(^estTgtPosit.y  -  PlatCenter.y),2.0)); 

DestBigZ  :=  DestTgtPosit.z  -  PlatCenter.z; 

END  IF; 

RangeCheck  :=  SQRT(POWER(curBigXY,2.0)  +  POWER(curBigZ,2  0)) 

REPEAT 

fiiture.x  :=  curTgtPosit.x  +  (SpeedVec.x*TimeInt);(New  x  pos} 
foture.y  :=  curTgtPosit.y  +  (SpeedVec.y*TinieInt);{Newy  pos} 
fiiture.z  :=  airTgtPosit.z  +  (SpeedVec.z*TimeInt);{New  z  pos) 
nextRange  :=  SQRT(POWER(future.x  -  PlatCenter.x,2.0)  + 

POWER(future.y  -  PlatCenter.y,2.0)  + 

PO\VER(future.z  -  PlatCenter.z,2.0)); 
curTgtPosh  :=  future; 

{  ***  Check  if  the  range  is  decreasing,  if  so,  add  to  array  if  R  <  Ro} 

IF  nextRange  <=  RangeCheck 

curBigXY  :=  SQRT(POWER((curTgtPosit.x  -  PlatCenter.x),2.0)  + 
POWER((curTgtPosit.y  -  PIatCenter.y),2.0)); 
curBigZ  :=  curTgtPosit.z  -  PlatCenter.z; 

RangeCheck  ;=  nextRange; 

IF  nextRange  <=  Ro 
counter  :=  counter  +  1 ; 

LookRanges[counter]  :=  nextRange; 

IF  counter  o  1 

TimeToRange[counter]  :=  (SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  POWER(future.y  -  RoPosit.y,2.0) 

+  POWER{future.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 
PO\VER(SpeedVec.y,2.0)  + 
POWER{SpeedVec.z,2.0))) 

+  TimeToRange[counter  -  1] 

RoPosit.x  :=  fiiture.x; 

RoPosit.y  :=  foture.y; 

RoPosit.z  :=  foture.z; 

ELSE 

{  ***  Time  to  go  from  Ro  to  first  accepted  position} 

TimeToRange[counter]  ;=  SQRT(POWER(foture.x  -  RoPosit.x,2.0) 

+  POWER(foture.y  -  RoPosit.y,2.0) 

+  POWER(foture.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 
POWER(SpeedVec.y,2.0)  + 
POWER(SpeedVec.z,2.0)); 

RoPosit.x  :=  foture.x; 

RoPosit.y  :=  foture.y; 

RoPosit.z  :=  foture.z; 

END  IF; 
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ENDEF; 

ELSE 

Inbound  :=  FALSE;{Tgt  passed  ovhd} 
curmanuv  :=  nextmanuv; 

END  IF; 

UNTIL  (nextRange  <  SQRT(POWER(DestBigXY,2.0)  +  POWER(DestBigZ,2.0)))  OR  (NOT(Inbound)); 
END  WHILE; 

{Backup  Plan  for  a  fly  through  situation} 

IF  Inbound 

RangeCheck  :=  SQRT(POWER(curBigXY,2.0)  +  POWER(curBigZ,2.0)); 

REPEAT 

fiiture.x  :=  curTgtPosit.x  +  (SpeedVec.x*TinieInt);{New  x  posit} 
future.y  :=  curTgtPosit.y  +  (SpeedVec.y*TimeInt);{New  y  posit} 
future.z  :=  curTgtPosit.z  +  (SpeedVec,z*TimeInt);{New  z  posit} 
nextRange  :=  SQRT(POWER(future.x  -  PIatCenter.x,2.0)  + 

POWER(iuture.y  -  PlatCenter.y,2.0)  + 

POWER(future.z  -  PlatCenter.z,2.0)); 
curTgtPosit  :=  future; 

{  Check  if  the  range  is  decreasing,  if  so,  add  to  anray  if  R  <  Ro} 

IF  nextRange  <=  RangeCheck 

curBigXY  :=  SQRT(POWER((curTgtPosit.x  -  PlatCenter.x),2.0)  + 

POWER({curTgtPosit.y  -  PlatCenter.y),2.0)); 
curBigZ  :=  curTgtPosit.z  -  PlatCenter.z; 

RangeCheck  :=  nextRange; 

IF  nextRange  <=  Ro 
counter  :=  counter  +  1 ; 

LookRanges[counter]  ;=  nextRange; 

{  ***  Time  to  go  fi’om  Ro  to  all  other  accepted  positions  } 

IF  counter  o  1 

TimeToRange[counter3  :=  (SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  POWER(future.y  -  RoPosit.y,2.0)  + 
POWER(fiiture.z  -  RoPosit.z,2.0)) 

/  (SQRT(POWER(SpeedVec.x  ,2.0)  + 
POWER(SpeedVec.y,2.0)  + 
POWER{SpeedVec.z,2.0))))  + 

TimeToRange[counter  -  1] ; 

RoPosit.x  ;=  fiiture.x; 

RoPosit.y  :=  future.y; 

RoPosit.z  ;=  future.z; 

ELSE 

{  ***  Time  to  go  fi'om  Ro  to  first  accepted  position} 

TimeToRange[counter]  :=  SQRT(POWER(future.x  -  RoPosit.x,2.0) 

+  POWER(future.y  -  RoPosit.y,2.0) 

+  POWER(future.z  -  RoPosit.z,2.0)) 

/  SQRT(POWER(SpeedVec.x  ,2.0)  + 
POWER(SpeedVec.y,2.0)  + 
POWER(SpeedVec.z,2.0)); 

RoPosit.x  :=  fiiture.x; 

RoPosit.y  ;=  future.y; 

RoPosit.z  :=  future.z; 

ENDEF; 

END  IF; 

ELSE 

Inbound  :=  FALSE; {Tgt  passed  ovhd} 
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curmanuv  :=  nextmanuv; 

END  IF; 

UNTIL  NOT(Inbound); 

END  IF; 

ELSE  {  ***  Check  the  next  set  of  points  (main  if  loop)} 
curmanuv  :=  nextmanuv; 
curTgtPosit.x  :=DestTgtPosit.x; 
curTgtPosit.y  :=  DestTgtPosit.y; 
curTgtPo^.z  :=DestTgtPosit.z; 

IF  curmanuv  o  NILOBJ 

nextmanuv  :==  ASK  TgtManuv  Next(curmanuv); 

END  IF; 

END  IF;  {main  if  loop} 

IF  NOT(Inbound) 

{  ***  The  array  is  stored  in  reverse  order  here} 

K  :=  counter;  (Lets  me  know  the  size  of  the  array} 

EXIT; 

END  IF; 

END  WHILE; 

DISPOSE(future); 

DISPOSE(SpeedVec); 

DISPOSE(TgtSpeed); 

DISPOSE{RoPosit); 

DISPOSE(DestTgtPosit); 

DISPOSE(PlatSpeed); 

DISPOSE(IPTgt); 

END  METHOD; 

ASK  METHOD  FindARange(IN  N  :  INTEGER) ;  REAL- 
VAR 

a,  b,  answer  :  REAL; 
i  :  INTEGER; 

BEGIN 

NEW(Ufour,  1..K); 

a  :=  (Threshold  -  FLOAT(N)  +  1.0)/FLOAT(N); 

{  ***  LookRmges  goes  from  [1]  =  furthest  out  range  from  platform[k]  =  closest  range  to  platform. 

The  ratio  of  R/Ro,  K  intervals.  This  is  an  array  of  U''4th  power} 

FORi;=lTOK 

Ufour[i]  :=  a  *  POWER((LookRanges[i]  /  Ro),  4.0); 

END  FOR; 

answer  :=  ASK  SELF  TO  FindSingleHitPd(N); 

DISPOSE(Ufour); 

DISPOSE(Pd); 

DISPOSE(CumProb); 

DISPOSE(Ratio); 

DISPOSE(LookRanges); 

DISPOSE(TimeToRange); 

RETURN(answer); 

END  METHOD; 

ASK  METHOD  FindSingleHitPd(IN  N ;  INTEGER)  :  REAL 
VAR 

i  ;  INTEGER; 
answer  :  REAL; 
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BEGIN 

NEW(Pd,  1..K); 

(Array  of  The  Single  Hit  Probability  ofDetection) 

FORi:=lTOK 

Pd[i]  :=EXP(-Ufour[i]); 

END  FOR; 

answer  :=  ASK  SELF  TO  BuildCDF(N); 

RETURN(answer); 

END  METHOD; 

ASK  METHOD  Bui]dCDF(IN  N  :  INTEGER) ;  REAL; 

VAR 

Check,  CumValue, 

a,  b,  CumSum,  c, 

hi, sum  :  REAL; 

answer  :  REAL; 

cdfrec  :  CDFRec; 

i,j  :  INTEGER; 

BEGIN 

NEW(Ratio,  1..K); 

{  ***  Find  the  R^o  ratio  and  store  it} 

FORi:=l  TOK 

Ratiop]  ;=  LookRanges[i]  /  Ro; 

END  FOR; 

{  ***  Now  Find  1  -  SHPD  and  load  to  matrix  } 

NEW(CumProb,  1..K); 

FORi:=lTOK 

NEW(CumProb[i]); 

CumProb[i].Range  :=  LookRanges[i]; 

{  ***  Find  the  Appropriate  Single  Hit  Pd  and  find  prob  no  detect  at  look  i) 
CumProb[i].NoDetect  :=  1.0  -  Pd[i];  (Probability  of  no  detection} 

END  FOR; 

(  ***  Finally  compute  the  cumulative  Pd} 

(  ***  Initialize} 

CumProb[l].CumPd  :=Pd[l]; 

CumSum  :=  CuniProb[l].CumPd; 

OUTPUTC'CDF  Range  Time  Index"); 

OUTPUT(CumSum,"  ",CumProb[l].Range,"  ",TimeToRange[l],"  1");  } 

FORi:=2TOK 
CumValue  :=  1.0; 

FOR  j  :=  1  TO  i 
IF  j  =  i 

CumValue  :=  CumValue  *  Pd[j]; 

ELSE 

CumValue  :=  CumValue  *  CumProb[j].NoDetect; 

END  IF; 

END  FOR; 

CumSum  :=  CumSum  +  CumValue; 

CumProbpJ.CumPd  :=  CumSum;  (Cummulative  Cdf  IS  BUILT} 

OUTPUT(CumProb[i].CumPd,"  ",CumProb[i] .Range,"  ",TimeToRange[i],  "  ",i) 

END  FOR; 

(  ***  In  case  CDF  does  not  go  to  1,  then  rescale  such  that  it  will  be  a  CDF} 

IF  (CumProb[K].CumPd  <  1.0  +  e)  AND  (CumProb[K].CumPd  >  1.0  -  e) 

FORi:=lTOK 
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CumProb[i].CumPd  :=  CumProb[i].CumPd  /  CuinProb[K].CumPd- 
ENDFOR; 

ENDEF; 

answer  ;=  ASK  SELF  TO  FindRandRange; 

RETURN(answer); 

END  METHOD; 

ASK  METHOD  FindRandRange ;  REAL; 

{^rse  transform  method  pg.  470  for  discrete  emperical  distribution} 

SetU,  P  :  REAL; 

I,i,j  ;  INTEGER; 

SelectedRange  ;  REAL; 

Foundl  :  BOOLEAN; 

BEGIN 

Foundl  ;=  FALSE; 

i  ;=  1; 

SetU  :=  ASK  RandVar  UniformReal(0.0,].0) 

REPEAT 

IF  (i  =  1)  AND  (SetU  <=  CumProb[i].CumPd) 

SelectedRange  ;=  TimeToRangep]; 

Foundl  :=  TRUE; 

ELSff  (SetU  >  CumProb[i].CumPd)  AND  (SetU  <=  CumProb[i+l].CumPd) 

SelectedRange  :=  TimeToRange[i+l]; 

Foundl  :=  TRUE; 

ELSE 

i:=i+l; 

END  IF; 

UNTIL  Foundl  =  TRUE; 

RETURN(SelectedRange); 

END  METHOD; 

END  OBJECT; 

END  MODULE. 

The  module  lAARONTIME.mod  stores  the  random  time  of  detection,  the  platform 
and  target  combination,  and  the  radar  sensor  system  that  is  detecting  the  target. 

Definition  module  DAARONTIME.mod; 

DEFINITION  MODULE  AARONTIME; 

FROM  RandMod  IMPORT  RandomObj; 

FROM  ListMod  IMPORT  QueueList; 

{  THE  PURPOSE  OF  THIS  MODULE  IS  TO  HOLD  ALL  ENTRY  DETECTION  TIMES} 

TYPE 

TimeRec  =  RECORD 


Plat 

:  STRING; 

Sensor 

.  STRING; 

Threat 

:  STRING; 

Time 

:REAL; 

END  RECORD; 
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TimeObj  =  OBJECT(QueueList) 
temp :  TimeRec; 

ASK  METHOD  Objlnit; 

ASK  METHOD  AddRecord(IN  P :  STRING;IN  S  ;  STRING; 

IN  T  ;  STRING;IN  Time  ;  REAL); 

ASK  METHOD  EmptyMe; 

END  OBJECT; 

VAR 

MasterEntryTime ;  TimeObj; 

END  MODULE. 


Implementation  module  lAARONTIME.mod: 

IMPLEMENTATION  MODULE  AARONTIME; 

FROM  RandMod  IMPORT  RandomObj; 

OBJECT  TimeObj; 

ASK  METHOD  Objlnit; 

BEGIN 

END  METHOD; 

ASK  METHOD  AddRecord(IN  P  ;  STRING;IN  S  ;  STRING;IN  T  ;  STRING;IN  Time  REAL) 
BEGIN 

NEW(temp); 
temp.PIat  :=  P; 
temp.Sensor  ;=  S; 
temp.Threat  :=  T; 
temp.Time  :=  Time; 

ASK  SELF  TO  Add(temp); 

END  METHOD; 

ASK  METHOD  EmptyMe; 

VAR 

ThisRec,  NextRec ;  TimeRec; 

I,  J  :  INTEGER; 

BEGIN 

J  :=  ASK  SELF  numberin; 

ThisRec  ;=  ASK  SELF  First; 

NextRec  :=  ASK  SELF  Next(lIiisRec); 

FOR  I  :=  1  TO  J 
IFIOJ 

ASK  SELF  TO  ReraoveThis(ThisRec); 

ThisRec  :=  NextRec; 

NextRec  :=  ASK  SELF  Next(ThisRec); 

ELSE 

ASK  SELF  TO  RemoveThis(ThisRec); 

END  IF; 

END  FOR; 

END  METHOD; 

END  OBJECT; 

END  MODULE. 
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The  module  IDetect.mod  is  the  engine  for  all  detections  that  take  place  in  a  given 
scenario.  It  uses  the  cookie  cutter  model  to  determine  the  time  of  detection.  Changes 
were  made  to  the  module  so  that  if  a  sensor  of  type  ‘RADAR’  is  being  evaluated  then 
apply  the  random  method  of  determining  a  detection  time.  If  not,  use  the  standard  cookie 
cutter  model.  The  listed  code  is  an  example  of  how  the  switching  between  random  and 
deterministic  routines  is  accomplished.  Similar  changes  were  also  made  in  the  method 
SolveCurvedProblem  of  IDetect.mod. 

Example  from  IDetect.mod  ASK  METHOD  StraightLineIntercept 

{****  Below  While  loop  by  ASE  15FEB96  it  performs  a  search  to  see  if  the 
detection  time  is  already  in  the  MasterEntiyTime  array.  If  it  is,  then  you  do 
not  want  to  recalculate  the  detection  time  if  not  needed  due  to  changes  in 
the  scenario.  ****} 

ThisRec  :=  ASK  MasterEntryTime  First; 

WHILE  (ThisRec  o  NILREC) 

IF  (ThisRec.Plat  =  ASK  Platform  ID)  AND  (ThisRec. Threat  =  ASK  Target  ID) 

{  Means  we  have  a  match  } 

IF  ThisRec.  Sensor  =  Sensomame 
RandTime  :=  ThisRec.Time; 

PrevAdded  :=  TRUE; 

OUTPUT("Have  a  match"); 

END  IF; 

ThisRec  :=  ASK  MasterEntryTime  Next(ThisRec); 

ELSE 

OUTPUT("  *♦*  Checking  MasterEntryTime  Array  No  match"); 

ThisRec  ;=  ASK  MasterEntryTime  Next(ThisRec); 

END  IF; 

END  WHILE; 

CASE  Solutions 

WHEN0;tl  := -1000.0; 

t2  := -1000.0; 

EntryExit.Entry.Time  :=  -1000.0; 

EntryExit.Exit.Time  := -1000.0; 

WHEN  1 :  tl  :=  ABS((TgtLWPoint.Location.x  -  Soil)  /  DeltaVel)  +  SimTime- 
t2:=tl; 

ASK  Environment  TO  SetEntryLocation(ASK  Target  PredictedPosition(tl)); 
sensorrec  :=  ASK  MasterSensorList  TO  GiveSensor(Sensomame); 

IF  (sensorrec.TypeSensor)  =  "RADAR" 

OUTPUT("Detected  a  RADAR  and  executing  random  routine"); 

ASK  Environment  TO  SetCurrentPlats(Platform,Target); 

IF  PrevAdded  =  FALSE  (Don't  want  another  rand  time} 

ASK  Environment  TO  FindDetectTime(sensorrec.Attributes.PulseInt,  ASK  Platform  ID, 
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sensorrec.  SensorTypeName, 

Target,  DetectionRadius); 

END  IF; 

DoRandom  :=  TRUE; 

END  IF; 

{If  random  routine  is  executed  do  this} 

IF  DoRandom 

IF  Prev  Added  =  FALSE 

EntryExit.Entiy.Time  ;=  tl  +  EnvironmentRandTime; 

DISPOSE(EntryExit.Entry.Location); 

OUTPUT("NEG  TI  found  in  Detect-straightKne  motion"); 

EntryExit.Entry .Location  :=  ASK  Target  PredictedPosition(tl  +Environment.RandTime); 
EntryExit.Exit.Time  :=  t2 ; 

DISPOSE(EntryExit.Exit.Location); 

EntiyExit.Exit.Location  :=  ASK  Target  PredictedPosition(tl  +Enviromnent.RandTime); 
OUTPUT("  Random  additional  time  to  intercept  in  IDetect: ", Environment. RandTime);’ 
ELSE 

EntryExit.Entiy.Time  :=  tl  +  RandTime; 

DISPOSE(EntryExit.Entry.Location); 

OUTPUT("NEG  Tl  found  in  Detect-straightline  motion"); 

EntryExit.Entiy.Location  :=  ASK  Target  PredictedPosition(tl+ RandTime); 
EntryExit.Exit.Time .- 12 ; 

DISPOSE(EntryExit.Exit.Location); 

EntryExit.Exit.Location  ;=  ASK  Target  PredictedPosition(tl  +RandTime); 

OUTPUT("  Random  time  ALREADY  IN  ARRAY  in  IDetect;  ",RandTimey 
END  IF; 

ELSE 

EntryExit.Entry.Time  ;=tl; 

DISPOSE(EntryExit.Entry.Location); 

OUTPUT("NEG  Tl  found  in  Detect-straightline  motion"); 

EntryExitEntry  Location  :=  ASK  Target  PredictedPosition(tl); 

EntiyExit.Ejdt.Time  :=  t2 ; 

DISPOSE(EntryEjdt.Exit.Location); 

EntryExit.Exit.Location  :=  ASK  Target  PredictedPositionltl) 

END  IF; 

WHEN  2: 

IF  (ABS(Distanc^TgtLWPoint.Location,PltLWPoint.Location))  - 
DetectionRadius  <  0.0  ) 
tl  :=  SimTime(); 

IF  DeltaVel  >  0.0 

t2  :=  ABS((TgtLWPoint.Location.x  -  MAXOF(Soll,Sol2))  /DeltaVel)  +  SimTime(); 

ELSE 

t2  :=  ABS{(TgtLWPoint.Location.x  -  MINOF(Soll,Sol2))  /DeltaVel)  +  SimTime()- 
ENDIF; 

ELSE 

IF  DeltaVel  >  0.0 

tl  :=  ABS((TgtLWPoint.Location.x  -MINOF(Soll,Sol2))/DeltaVel)  +  SimTirneQ; 
t2  :=  ABS((TgtLWPoint.Location.x  -MAXOF(Soll,Sol2))/DeltaVel)  +  SimTimeO; 

ELSE 

tl  ;=  ABS((TgtLWPoint.Location.x  -MAXOF{Soll,Sol2))/DeItaVel)  +  SimTirneQ; 
t2  :=  ABS((TgtLWPoint.Location.x  -MINOF(Soll,Sol2))/DeltaVel)  +  SimTirneQ" 

END  IF; 

IF  t2  <  tl 
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t2  :=tl; 

END  IF; 

END  IF; 

ASK  Environment  TO  SetEritiyLocation(ASK  Target  PredictedPosition(tI)); 
sensorrec  ;=  ASK  MasterSensorList  TO  GiveSensor(Sensomame); 

IF  (sensorrec.TypeSensor)  =  "RADAR" 

OUTPUT("Detected  a  ",  sensorrec.TypeSensor, "  which  should  match  RADAR  and 

.  .  executing  random  routine"); 

ASK  Environment  TO  SetCuiTentPlats(Platfonn, Target); 

IF  PrevAdded  =  FALSE  {Don’t  want  another  rand  time) 

ASK  Environment  TO  FindDetectTime(sensorrec.Attributes.PulseInt,  ASK  Platform  ID, 

sensorrec.  SensorT3rpeName, 

Target,  DetectionRadius); 

END  IF; 

DoRandom  :=  TRUE; 

END  IF; 

{If  random  routine  is  executed  do  this } 
ff  DoRandom 

IF  PrevAdded  =  FALSE 

EntryExit.Entry.Time  :=  tl  +  Environment.RandTime; 
DISPOSE(EntryExit.Entry.Location); 

EntiyExit.Entry  Location  :=  ASK  Target  PredictedPosition(tl  +Environment.RandTime)- 
DI SPO  SE(EntryExit .  Exit .  Location); 

EntiyExit.Exit.Time  :=  t2; 

EntryExit.Exit.Location  :=  ASK  Target  PredictedPosition(t2); 

OUTPUTC  Random  additional  time  to  intercept  in  IDetect:  ".Environment.RandTime); 
ELSE 

EntryExit.Entry.Time  ;=tl  +RandTime; 

DISPOSE(EntryExit.Entry.Location); 

OUTPUT("NEG  Tl  found  in  Detect-straightline  motion"); 

EntryExit.Entiy.Location  ;=  ASK  Target  PredictedPosition(tl+  RandTime); 
EntiyExit.Exit.Time  :=  t2 ; 

DISPOSE(EntryExit.Exit.Location); 

EntryExit.Exit.Location  :=  ASK  Target  PredictedPosition(tl  +RandTime); 

OUTPUTC  Random  time  ALREADY  IN  ARRAY  in  IDetect:  ".RandTime) 

END  IF; 

ELSE 

EntryExit.Entry.Time  :=  tl; 

DISPOSE(EntryExit.Entry.Location); 

EntryExit  Entry-Location  ;=  ASK  Target  PredictedPosition(tl); 
DISPOSE(EntryExit.Exit.Location); 

EntryExit.Exit.Time  :=  t2; 

EntryExit.Exit.Location  :=  ASK  Target  PredictedPosition(t2) 

END  IF; 

OTHERWISE 

WriteLine{"Fuimy  return  from  Quadratic"); 

END  CASE; 

OUTPUTC  Tl  intercept;  ",tl); 

OUTPUT("  The  time  ",ASK  Platform  ID,"  will  detect  ",ASK  Target  ID,  "  =  ".EntryExit.Entry.Time); 
{Want  to  take  a  data  point  of  the  random  time  until  detection.  This 
stat  could  also  be  gathered  on  only  the  random  times  generated  earlier. 

Gathering  stat  here  combine?  any  changes  in  the  initial  time  Ro  is 
reached ...} 

StatTime  ;=  EntryExit.Entry.Time; 
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ASK  StatKeeper  TO  GetSaniple(StatTiine); 

OUTPUT("Gathered  a  Stat"); 

END  METHOD; 

Many  other  changes  were  made  to  the  NPS  Platform  Foundation  that  will  not  be 
discussed  in  this  appendix.  The  changes  insured  the  flow  of  the  correct  information  to 
the  proper  modules.  This  information  exchange,  althou^  initially  done  only  for  this 
thesis,  allows  for  many  other  uses  of  the  NPS  Platform  Foundation.  All  new  and 
modified  modules  for  this  thesis  are  listed  below.  An  asterisk  denotes  new  modules. 

•  D/IAARONLmod  * 

•  D/IAARONCDF.mod  * 

•  D/lAARONTIME.mod  * 

•  D/ICPAmod 

•  D/lDetect.rriod 

•  D/IDetRng.mod 

•  D/IMathFunc.mod  * 

•  D/IPList.mod 

•  D/lReact.mod 

•  D/ISpace.mod 

•  D/IVS.mod 

•  D/IVSMan.mod 

•  D/IWrite.mod 

•  D/IGP.mod 

•  IGPMast.mod 
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•  IGSimExc.mod 


•  IMamcv.mod 

•  IRep.mod 
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